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ABSTRACT 

We  construct  an  Isotropic  approximation  to  the  Pokker-Planck 
collision  operator  (IPP  model)  for  the  purpose  of  studying  the 
small  amplitude  wave  motion  of  an  unbounded,  fully  ionized,  neutral, 
two  component  plasma  in  a  steady  applied  magnetic  field.   We  solve 
the  linearized  Boltzmann's  equations  with  Maxwell's  equations  and 
construct  the  dispersion  relation.   Because  the  isotropic  Pokker- 
Planck  operators  describe  a  diffusion  process  in  velocity  space 
while  the  previously  considered  Krook  relaxation  models  do  not 
we  find  an  infinite  number  of  new  Larmor  resonances  modes  at 
arbitrary  directions  of  propagation  with  respect  to  the  applied 
magnetic  field.   The  new  resonances  remain  in  the  limit  of  long 
macroscopic  wavelengths  and  small  collision  frequencies.   At 
propagation  parallel  to  the  applied  magnetic  field  we  find  new 
resonances  for  wavelength  large  and  small  with  respect  to  the  mean 
free  path. 

In  other  macroscopic  ranges  the  Krook  and  IPP  models  give 
similar  results.   As  a  special  example  we  verify  the  Paraday 
rotation  formula  for  a  particular  collisional  plasma 
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1.   Introduction 

We  Introduce  an  Isotropic  approximation  to  the  Pokker-Planck 
collision  operators  In  Boltzmann's  equations  In  order  to  study 
small  amplitude  wave  motion  In  a  fully  Ionized,  neutral,  two 
component  plasma  In  a  steady,  applied  magnetic  field.   With  this 
model  we  find  Larmor  resonance  modes  which  can  not  be  found  from 
theories  based  on  the  relaxation  model  of  Bhatnagar,  Gross  and 
Krook    (the  Krook  model).   We  remark  that  so  far  the  Krook  model 
Is  the  only  other  collision  model  from  which  one  computes  Larmor 
resonances.   Moreover,  these  modes  persist  In  a  certain  non-uniform 
colllsionless  limit. 

The  complexity  of  the  collision  operators  either  In  Boltzmann 
Integral  form  or  In  Pokker-Planck  form  forces  one  to  seek  simpler 

models  to  study  wave  motion  In  plasmas  from  a  microscopic  approach. 

[pl 
The  model  of  R.L.  Llboff    Is  a  generalization  of  the  Krook  model 

which  Includes  the  effects  of  Ion-electron  collisions.   His  model, 

therefore,  has  the  attractive  feature  of  giving  the  macroscopic 

conservation  equations  when  one  computes  the  required  moments  of 

Boltzmann's  equations.   A  defect  of  all  Krook  models  Is  that  they 

do  not  describe  a  diffusion  process  In  velocity  space.   The 

Isotropic   Pokker-Planck  (IPP)  model  Introduced  here  has  two 

Important  features.   It  gives  the  macroscopic  conservation 

equations  and  It  describes  a  diffusion  process  In  velocity  space. 

Both  properties  are  shared  with  the  Boltzmann  Integral  operators 

[3l 
and  the  original  Pokker-Planck  operators    .   Por  those  effects 
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which  are  sensitive  to  the  collision  model  we  expect  to  find 
answers  different  from  those  found  from  Krook  model  theories.   We 
verify  that  the  behavior  of  Larmor  resonances  found  by  the  IFF 
model  is  radically  different  from  the  behavior  of  those  found  by 
Llboff    .   On  the  other  hand  In  certain  parameter  ranges  there  are 
effects  which  are  not  sensitive  to  the  collision  model.   An  example 
is  given  by  the  plasma  electromagnetic  and  magnetohydrodynamic 
modes  in  the  long  wavelength  limit.   Liboff  shows  that  to  lowest 
order  the  microscopic  and  macroscopic  description  of  these  modes  are 
equivalent.   The  IPP  results  agree  to  lowest  order  with  Llboff 's. 

■X- 

We  construct  the  dispersion  relation  from  the  formal  solution 
to  the  linearized  Boltzmann's  equations  and  Maxwell's  equations. 
The  roots  of  the  dispersion  formula  are  connections  between 
frequency  and  waveniimber  and  are  Interpreted  variously  as 
oscillations,  resonances,  propagating  modes  and  the  like.   In  the 
detailed  analysis  of  the  roots  we  prescribe  that  the  wavenumber 
vector  be  real  and  that  the  wavelength  be  large  or  small  with 
respect  to  certain  representative  physical  lengths.   We  then 
compute  frequency  in  terms  of  wavelength  and  the  physical  parameters 
from  the  dispersion  equation. 

Our  principal  results  are  in  the  study  of  the  Larmor 
resonances.   These  resonances  are  roots  of  the  dispersion  equation 
for  which  the  frequency  is,  roughly,  an  Integral  multiple  of  the 


The  terms  "Hisperson  equation",  "dispersion  formula"  and 
"dispersion  relation"  are  used  Interchangeably. 
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Larmor  frequency.   There  is  a  radical  difference  in  the  kind  of 

resonances  predicted  from  the  IPP  dispersion  equation  and  from 

rol 
those  predicted  from  the  dispersion  equations  of  the  Krook  model 

r  -1  o  "I 
or  from  the  colllslonless  Boltzmann's  equation     . 

Prom  the  Krook  and  colllslonless  Boltzmann's  equations  one 

finds  an  Infinite  number  of  Larmor  resonances  only  at  propagation 

exactly  perpendicular  to  the  applied  magnetic  field.   Prom  the  IFF 

equation  one  finds  an  Infinite  number  of  Larmor  resonances  at 

arbitrary  angles  of  propagation  (except  parallel)  to  the  applied 

magnetic  field.   This  reflects  the  physical  property  that  the  IPP 

model  describes  a  diffusion  process  In  velocity  space  while  the 

Krook  model  and  the  colllslonless  Boltzmann's  equation  do  not. 

There  Is  another  difference  between  the  Krook  and  IPP  models.   At 

perpendicular  propagation  ,'  to  each  Krook  resonance  there  Is  an 

Infinite  number  of  IPP  resonances  separated  by  Integral  multiples  of 

the  collision  frequency.   This  reflects  the  mathematical  property 

that  the  IPP  differential  operator  Is  unbounded  with  a  discrete 

operafoV' 

spectrum  while   the   Krook   flqiii^t'TQjT  Is  bounded. 

The  distinctions  between  the  two  colllslonal  models  remain  In 
the  non-uniform  colllslonless  limit  obtained  by  keeping  collision 
length  to  wavelength  fixed  and  small  as  the  collision  frequency 
becomes  arbitrarily  small. 

A  precursor  to  the  IPP  operator  Is  the  Pokker-Planck  operator 

[41  [Sl 

familiar  from  Brownlan  motion  studies    .   A.  A.  Logunov  -^  used 

such  an  operator  to  study  the  space  periodicity  of  gas  discharges. 
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[61 
Somewhat  later  A.  Lenard  and  I.  Bernstein    used  the  same 

operator  to  study  plasma  oscillations.   Neither  of  these  authors 

considered  the  problem  for  a  steady  applied  magnetic  field. 

I  7    o  "I 

However  J.  Keilson,  S.  Schneider  and  A.  Cantor  '' '    did  study 
propagation  in  the  upper  atmosphere  in  the  presence  of  a  steady 

magnetic  field.   Their  results  are  contained  among  the  lowest  order 

[21 
roots  found  by  Liboff    .   In  none  of  these  models  were  the  ion- 
electron  interactions  included  so  that  the  macroscopic  equations 
of  motion  of  fluid  dynamics  are  lost. 

In  section  2  we  describe  the  construction  of  the  isotropic 
Fokker-Planck  model  from  the  full  Pokker-Planck  equation.   Following 
this  the  equations  of  motion,  Boltzmanri's  equations  and  Maxwell's 
equations,  are  linearized  about  an  absolute  equilibrium  state. 

Secton  3  contains  the  formal  solution  to  the  linearized 

equations  and  the  construction  of  the  dispersion  relation.   At  the 

[2l 
end  of  this  section  we  show  how  Libof f ' s    lowest  order  results 

may  be  obtained  from  our  equations. 

Section  4  contains  the  principal  results  of  this  work.   These 

are  the  studies  of  the  Larmor  resonances  in  different  limiting 

cases.   Following  the  introductory  section  4A  we  describe  in  4B 

the  analysis  followed  in  the  rest  of  the  section  except  in  the 

short  wavelength  limit.   Here  too  we  indicate  the  reason  why  both 

the  collisional  Krook  model  and  the  collisionless  model  give  an 

infinite  number  of  Larmor  resonances  only  at  angles  of  propagation 

exactly  normal  to  the  direction  of  the  applied  magnetic  field. 


Section  4C  contains  our  study  of  the  non-uniform  limit.   Section 
4D  contains  the  uniform  colllsionless  results  at  angles  of 
propagation  normal  to  the  applied  magnetic  field.   In  section  4E 
we  consider  the  collisional  Larmor  resonances  at  arbitrary  directions 
of  propagation,  at  propagation  normal  to  the  magnetic  field  to 
compare  with  Llbof f ' s  results,  and  at  propagation  parallel  to  the 
magnetic  field.   All  of  these  results  are  in  the  long  wavelength 
limit,  that  is,  long  with  respect  to  the  Larmor  radius  and  long 
with  respect  to  the  collision  mean  free  path  except  for  the  uniform 
colllsionless  case  where  the  concept  of  mean  free  path  does  not 
enter.   In  section  4P  we  present  our  results  for  propagation 
parallel  to  the  magnetic  field  for  wavelength  short  compared  with  a 
mean  free  path. 

In  section  5  we  first  calculate  the  effect  of  a  plasma  in  a 
magnetic  field  on  the  electromagnetic  propagation  mode  for 
propagation  parallel  to  the  applied  magnetic  field.   The  calculation 
is  given  for  a  set  of  experimental  parameters  suggested  by  Dr.  Fred 
L.  Rlbe.   In  doing  this  calculation  we  verify  the  classical  formula 
describing  the  Faraday  rotation  in  the  parameter  range  of  Dr. 
Ribe's  experiments. 

We  then  calculate  the  effect  of  collisions  on  the  plasma 
oscillation  mode  and  compare  the  results  with  those  of  previous 
authors. 

Section  6  contains  the  summary  and  conclusions  of  this  study. 
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2.   Equations  and  Linearization 
A.   Equations 

We  start  with  the  Boltzmann's  equation  In  which  the  collision 
processes  are  represented  by  Pokker-Planck  differential  operators  . 
These  operators   J  have  the  typical  form 

a2    w.    „,       a 


J(f)  =  i(^^^)Ks^)  -  gr  (V) 


r   s  r 


■X-* 


where   f   Is  the  single  particle  distribution  function,   ^    Is 

the  microscopic  particle  velocity  with . cartesian  components   ^^, 

b    and   a   are  the  diffusion  and  friction  coefficients 
rs        r 

respectively.   The  coefficients   b    and   a   are  In  turn 
complicated  Integrals  In  velocity  space  Involving  the  distribution 
function  of  the  scattering  particles. 

The  construction  (not  derivation)  of  the  Isotropic  Pokker- 
Planck  (IPP)  model  proceeds  as  follows.   Plrst  we  evaluate   b 
and  a   as  if  each  species  of  scattering  particles  were  In  a 
local  Maxwelllan  state.   Then  we  evaluate  the  coefficients  in  the 
limit  in  which  the  particle  energy  is  much  smaller  than  the 
average  energy  of  the  scatterers.   Finally  we  manipulate  the  coef- 
ficients involving  ion-electron  collisions  until  the  following 
characteristics  of  the  original  Pokker-Planck  coefficients  reappear, 
The  collision  terms  must  vanish  when  both  species  are  in  the  same 


*   c.f.  L.  Landau  [9],  K.Grad  [lO],  Rosenbluth  et  al.  [ll], 
**  The  wavy  underline  (,vw)  Indicates  the  vector  quantities. 
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local  Maxwellian  state.   The  conservation  equations  for  mass 
density,  total  momentum  density  and  total  energy  density  must  be 
satisfied  .   As  a  consequence  an  H-theorem  for  the  gas  mixture 
Is  also  satisfied.   We  now  postulate  that  the  coefficients  are 
valid  throughout  the  entire  velocity  space,  and  call  the  new  dif- 
ferential operator  the  Isotropic  Fokker-Planck  operator.   The 
results  are  given  In  equations  (2.6)  and  (2.?)  below. 

Before  writing  the  equations  we  need  to  Introduce  some 
notation  and  definitions.   The  gas  under  study  Is  a  fully  Ionized 
mixture  of  electrons  and  Ions  of  one  species.   Ion  quantities  carry 
a  caret  (^)  while  electron  quantities  remain  bare.   The  single 
particle  distribution  functions   f  and   f   are  functions  of  time   t, 
space  X,   and  velocity   |.   In  equations  (2.1)  we  define  the 
following  moments  of   f   and   f  :   The  number  densities  In 
configuration  space   n  and  n  ,  the  macroscopic  velocities   u  and 
u,   and  the  macroscopic  temperatures   T  and  T  . 

n  =  Jf6?i  ,  n  =  /fd^^  , 

nu  =  J^fd?  ,  nu  =  ftfd^l,  (2.l) 

3nRT  =  /(^-u)^fd^^      ,    3nRT  =  /(e-u)^fd^^. 


where  R  =  x/m,   R  =  K/m,   k  Is  Boltzmann's  constant  and  m  and  m 
are  the  particle  masses. 


*   Chapman  and  Cowling  [12]  Chapter  8. 

**  Chapman  and  Cowling  [12]  Chapter  4,  and  H.  Grad  [13] 
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In  equations  (2.2)  we  define  the  total  number  density  n,   the 

total  mass  density  p",   the  total  macroscopic  velocity  _u,   the 

I       '^  I  — 

temperatures   T   and  T   measured  with  respect  to  _u_,   and  the 

total  gas  temperature   T  . 

n=n+n,  p=mn+mn, 

/\  /\  /\ 

p  u  =  mnu  +  mnuj,  (2.2) 

n  i.l'-  =  nT  +  nT   , 

Finally  the  total  charge  density  q   and  current  density   j 
are  defined  In  terms  of  the  signed  electronic  and  Ionic  charges   e 
and  e   In  equations  (2.3)' 


q  =  ne  +  ne  , 


J  =  neu  +  neu 


(20) 


The  equations  of  motion  are  Maxwell's  equations  (2.4)  and  the 
two  Boltzmann's  equations  (2-5). 

V-E  =  q/e   ,     V-B  =  0  , 

S7^E=   -(a/at)^,  (2.4) 

VxB=  ^^1+  c-2(a/st)^, 

where  _E  and  _B_  are  the  electric  and  magnetic  field  vectors 
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respectively,   e   and   ll   are  the  dielectric  constant  and 
permeability  of  free  space  respectively,  and   c  =  (e  [x  )~^      Is  the 
speed  of  light.   We  use  rationalized  MK-§>  units  throughout. 


dt   'v?^  Y     ni'~-  ^  ^■'     ^  ?     Mt^coll 


m      ^- 


(2.5) 


where  \7e   is  the  gradient  operator  in  velocity  space,  and  the 

collision  terms   (xt")   n -,   and   (31-)   nn   are  now  those  of  the 

ot  coll       Mt'coll 

Isotropic  Fokker-Planck  model: 


(i)ooll   =   Jl(f)    -  J2(f). 


(2.6) 


/\  /\  /s  /\ 


(||)coii  =  Ji(f)  +  Jatf). 


where 


Jl{f)  =  V,  5|-(RT||-+(e^-u^)f), 
J,(f)=v^3|-(Rf|f-+  (5,  -u,)f), 
Jl(f )   =   v^  ^  (RT  If  +    (?^   -   u^)f). 


(2.7) 


/\  /S 


"p  I* 


The  collision  frequencies   v-,  ,  v-,   and   Vp   are  unknown  functions 
of  space  and  time.   That  they  are  unknown  is  the  penalty  exacted 
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for  simplifying  the  Pokker-Planck  coefficients.   The  collision 
terms  have  the  following  properties.   They  vanish  when   f   and   f 
are  In  the  same  local  Maxwelllan  state,  that  Is  when 


f   =  n(27rRT)      ^^^    exp(  -  [^-u]  VsRT)    , 


/\  /\  /\/\ 


f  =   n(27rRT)~   ^/^    exp(  -  [^-u]  ^/2RT)    , 


/Vv|    ^IVj 


u  =   u  and      T  =   T 


The  macroscopic  conservation  equations  follow  from  the 
Identities  (2.8). 


\J 


yj 


J^(f)d^|   =      J,(f)d^^   =    0, 


^1       /\      /v 


m    le    J^(f  )d^|   +  m   U    oT    (f  )d^^ 


=    0, 


|pj2(f)d^^   +1   L2j^(f)d5^   =   0. 


(2.8) 


In   addition   one   can   derive   an  H-theorem.      We   define 

H  =  /f   log  f  d^^   +  /f   log  f  d^|.      The   contributions   to   ^^ 


^     coming 
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from  the  collls^lon  terms   (37^)   -1-1  <  0,   with  equality  when   f   and 
f   are  In  the  same  local  Maxwelllan  state.   Finally  the   IPP 
operators  describe  a  diffusion  process  in  velocity  space,  a 

property  shared  with  the  Boltzmann  Integral  operators  and  the  full 

r-5] 
Pokker-Planck  operators   ^ 


B.   Linearization 

We  linearize  the  equations  of  motion  about  the  equilibrium 
solutions   f   and   f   in  the  rest  frame  of  the  total  gas.   This 
means  that  all  of  the  quantities:   the  distribution  functions,  their 
moments,  and  the  field  variables  are  written  in  the  form 


X  ^x^  +  X  ,  (2.9) 

where  the  arrow  signifies  that  wherever  X  appears  in  the  equations 

and  definitions  it  is  replaced  with  the  sum  of  its  equilibrium 

value  X   and  its  perturbation  value   X.   Products  of  two  or  more 
o 

perturbed  quantities  are  dropped  systematically.   The  equilibrium 
solutions  are  the  absolute  maxwelllans 

^o  =  "o^^ttC^)-  5/2  g^p(_^2/2c2)^ 

(2.10) 
f^  =  n^(27rc2)-  ^/^  exp{-^^/2C^ ) , 

where      C^   =   RT    ,      C^   -   RT      ,      T     =   T    ,      and   the   quantities      n    , 
o  000  o 

n        and     T        are    constants.      The   property  that   the   equilibrium 

^   /\ 
plasma,     be     neutral  means   that      fe[n     =    [e|n^. 
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We  suppose  for  simplicity  that  the  positive  particles  are 


singly  ionized  so  that  from  this  point  on  n 


o 


n 


Once  we  have 


linearized  there  is  no  difference  between   T  and   T   and  between 
T  and  T   since  they  differ  by  second  order  quantities  in  the 
perturbed  variables.   The  full  list  of  variables  appears  in 
equations  (2.II). 


f  -^  n^(g^  +  g) 

n  — ^  n   +  n 

^     o 


u  — >  0  +  u 


T  — >  T   +  T 

'      o 


q  -^  0  +  q 


E  ^  0  +  E 


where  ng  =f,ng  =f 
00    o"^   00    c 


f  -^  "o^^o  ^  ^^' 


P  ^  P^  +  P 


o 


u  — ^  0  +  u  , 

T  — ^  T  +  T, 
^   o     ' 


n  — ^  n   +  n  , 
-^  o 


u  — ^  0  +  u  , 


T  -^  T^  +  T  , 
1"  -^0  +  J  , 
B  -^B   +  B  , 


the  perturbed   f  =  n  g. 


the 


(2.11) 


perturbed   f  =  n  g,   and  _B 
^  o^       '^o 


is  the  constant  applied  magnetic 
field.   In  terms  of   g  the  moments  are  given  by 


n  =  n 


gd^^  , 


u  = 


(T/T  )  ^    \[i^/3C^    -    Ugd^^, 


(2.12) 


and  similarly  for  the  ion  quantities. 

Maxwell's  equations  (2.4)  remain  unchanged  in  form. 
Boltzmann's  equations  (2.5)  and  the  collision  terms  (2.6)  and 
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(2.7)   become 


If  +  4-Vs  +i,^£,-VeS  -  Jo(8)  =  l^+  ^  'S  -  5'i-o 


If  +i;V8  -ix&Ves  -  J„(g)  =  )^+  2(4  -  3)1  g„. 


C^  C 


where 


j„(e)  =  V  gf  (c^  I^  +  e^g)  , 

r  r 


J„(g)    =   V  ^   (02  1^  .   t^g)    , 

r  r 


e    T-,  — 

a  =  —  E  +   v-,u  +   v„u    , 

-m     --VV,  lv^-,   •  2'^. 


b 


(v^T   +    V2T)/T^      , 


/\  /\        /N 


b   =    (v^T   +    v^T)/T^      , 


/\  /\  /\ 


fi=eB/m      ,      n=eB/m, 


/\  /s 


(2.13) 


V    =     V-|_    +    Vg     ,        V    =     V-j^    +    Vg 


The  collision  frequencies   v-,  ,      v-,   and   Vp   are  unknown 
constants  In  this  theory.   None  the  less,  it  is  useful  to  have 
expressions  for  them  which  indicate  dependence  on  physical 
parameters  such  as  mass  and  temperature.   For  this  reason  we  give 
in  equations  (2.l4)  expressions  for   v-,   and   v-,  which  come 
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directly  from  the  original  Pokker-Planck  coefficients,  and  an 

expression  for   v„   which  derives  from  the  time  of  equlpartltlon 

of  energy  discussed,  for  example.  In  chapter  5  of  L.  Spltzer's 
[14]^ 

e   n    log/\ 
_  ___^ o 

1        3m2    £2(2TrKT    )5/2 


book 


o'  o 


^  =^  (f  ri  ' 


(2.14) 


m 


Smm^ 
"2  °  i;:^?72   ^1 


where 


A 


-(%)*  ff  °^''^ 


We  assume  that  equations  (2.14)  give  the  correct  dependence 
of  the  colllsilon  frequencies  on  the  Ion  and  electron  masses. 
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3 .   Construction  of  the  Dispersion  Relation 
A.   Program  of  the  construction 

This  brief  outline  of  the  various  steps  which  lead  from  the 
equations  of  motion  to  the  final  dispersion  relation  Is  Intended 
partly  as  a  guide  to  the  details  and  partly  as  a  way  of  defining 
what  we  mean  by  a  dispersion  relation  In  terms  of  Its  construction. 

Subsection  B  contains  the  formal  solution  to  the  Boltzmann's 
equations  (2.13)-   We  Illustrate  with  the  equation  for   g  which 
can  be  written 

L(t,x,  ^)g  =  R  g   . 

L  Is  the  differential  operator;   R   Is  a  linear  combination  of  E 
and  the  unknown  moments   u,  u,  T,   and   T.   Under  the  Laplace- 

r  -     ?  ■         r  •    ? 

Fourier  transformation    I  e   dt   e  -^'^   d  x   e  ■^'^d'^i        the 

'Oq       J  J 

equation  becomes 

L(s,k,w)g  =  R  g  +  g 

where   g.   Is  the  Initial  value  of  g,      and  L  Is  now  a  first 
order  linear  partial  differential  equation  In  w.   We  Invert  the 
equation  and  get 

g(s,k,w)  =  L"  (R  g  +  g. )  . 

This  means  that   g   Is  known  In  terms  of   E  and  the  moments  u, 
u,   T  and  T.   There  Is  a  similar  equation  for   g.   Maxwell's 
equations  and  the  equations  of  continuity  enable^  us  to  write 
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E  (s.k)   In  terms  of  the  moments   u(s,k)   and  u(s,k)   and  the 
initial  values  of  the  field  quantities.   The  equations  for   g  and 
g  are  not  closed.   That  is,   g  and  g  are  given  in  terms  of  their 
unknown  moments  u(s,k),   u(s,k),   T(s,k)   and   T(s,k),   We  cannot 
know  g  and   g  unless  we  know  those  moments.   In  subsection   C 
we  close  the  equations  and  solve  formally  for  the  moments.   This 
procedure  gives  us  the  dispersion  relation  which  is  our  goal. 

The  construction  given  in  subsection  C  rests  on  the  identities 

n(s,k)  -^  g(s,k,0)  , 


u(s,k)'NyV„g(s,k,w) 


W=0 


etc.   With  these  relations  and  the  corresponding  ones  for  the  ion 
variables  we  can  form  the  relation 

A(s,k)  U(s,k)  =  V(s,k) 

where  A   is  now  a  known  matrix,   U  is  a  column  vector  consisting 
of  E  and  the  moments   u,u,   T  and  T,   and  V  is  a  column  vector 
containing  the  initial  data.   The  dispersion  relation  is 


det  A(s,k)  =  0. 

The  roots  of  the  dispersion  relation  are  equations  connecting 
s   and  k.   It  is  clear  that  knowledge  of  the  analytic  properties 
of  the  initial  data  will  tell  us  which  of  the  roots  of  the  dis-  - 
persion  relations  are  relevant  to  a  given  problem^-^..   We  will 
neglect  this  aspect  of  the  problem  and  concentrate  on  the  dis- 
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perslon  relation  which  contains  all  possible  roots. 

B.   Formal  Solution  to  the  Linearized  Boltzmann's  Equations 
Under  the  transformation 

00 


O 


r  r 

J         ^ 


dt  I  e '  d  X  1  e  -v^-X  d^  (3-1 ) 


the  first  of  equations  (2.13)  becomes 


(s  +  vC  w  )g(s,k,w)  +  (vw  -  k  +  wxfi) -XZ  g( s,k,w) 


(1  a(s,k)-w  -  c2w2b(s,k))e"^  w  /2  ^  g.(k,w). 


(3.2) 


where   g. (k,w)   is  the  transform  of  the  initial  value  of   g  at 
time   t  =  0.   The  equation  for   g   is  similar  with  the  caret  (  ) 
over  the  appropriate  quantities.   Although  we  will  no  longer  carry 
the  Initial  values  in  the  transformed  Boltzmann's  and  Maxwell's 
equations  their  presence  is  Implied  throughout  the  development. 
To  integrate  equation  (3-2)  we  look  for  solutions  to 
d  w(t  ) 


'^    -  vw(t)  -  k  +  w(T))fn  (3.3) 


subject  to 


w(0)  =  w  .  (3-^) 


If  w(t)   satisfies  (3-3)  we  can  write  equation  (3.2)  in  the  f 


orm 
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,2  2, 


(s+vcV(T))g  +  ff  =  (la.w(T)  -  cV(T)b)e-C  ^  (^)/2  ^     (3^3) 


Equation  (5.5)  can  be  Integrated  easily.   The  condition  {^A) 
means  that  we  can  Integrate  equation  (3-5)  fo^  all  ^w. 
Solutions  to  (3.3)  are  of  the  form 


-T-^ 


w 


where  \i(-^\}  \^fn\'  ^'^'^  ^f^\  ^^'^  vectors-  Independent  of  x,  and 
a  Is  a  constant  to  be  determined.  We  Impose  condition  (3-^)  and 
get,  after  some  straightforward  vector  manipulation. 


w 


y  =  ji(i)  ^  -""' )  (ri(i))-£/"'  -  -''":^.  ^  -"'V  i       (3.6) 


-A 


where 


.2,  2 


^i/,N  =  [kv   +  vftyk  +  fi(f2-k)]/v(fi  +v  )   , 

fi  =  e  B  /m  carries  the  sign  of  the  electron  charge  e. 


and 


7 


+ 


fi 


I 


^2      -  1 


The  Integrating  factor  for  equation  (3-5)  is 


exp 


JT, 


[s+vC  w  (t  )]dT  .   The  lower  limit   t    Is  arbitrary  and 


drops  out  In  the  subsequent  manipulations.   We  multiply  equation 
(3-5)  by  the  Integrating  factor  and  Integrate  from  t  =  -00  to 
T  =  0.   The  lower  limit  on  the  left  hand  side  vanishes  In  the 
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range  of  validity  of  the  original  transformation  (3.1).   The 
result,  after  some  trivial  changes  in  the  variables  of  integration, 
is 


00 


g(s,k,w)  = 


\J 


dy  e^(y'^)  H(y,w), 


(3.7) 


where  j6{y ,v)      and  H(y,w)  '  are  defined  In  equations  (3-8)  and 


(3-9)  respectively. 
j2^(y,w)  = 


(s+vC   ^i^))y   -  — 


p2    2 
C   w 


C      (w-^L/       N  ) 


where 


-  r  e-^^(i-e-^^-^'^^y)  , 


(3.8) 


+liA 


(v+ifi)^ 
v^+f2^ 


and 


1 

2 


fi^ 


H(y,w)   =   1  U^^)-a^+    (^-b^l))- 


(fi-a) 
f2  l:n_:V   e-^y  +  A^  e- 


'+if2)y 


a' 


+  A      e 


■(v-ifi)y 


-w2(-y)c2b    , 


where     A.     =  -^ 


>kv^  I 


(3.9) 


fi 


~        1 


n 
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w 


^(-y)  =  ^if^^  +  2(w-y-(i)) 


2 


'■^   -^    A^{   X   I 


)^ 


-vy 


Q.' 


p^-(v+in)y 


+  I  ,e 


and 


r< 


+   e 


-(v-lfi)y 


k-ii,(l)^  ^ 


2  -2vy 


M. 


.2  2 


2\  /   2/  2,^2 


Equation  (3.?)  Is  a  formal  solution  to  Boltzmann's  equation. 
That  Is,   g(s,k,w)   Is  given  linearly  In  terms  of  the  unknown 

/\  /s  ^ 

moments   n(s,k),  n(s,k),  u(s,k),  u(s,k),  T(s,k)  and  T(s,k),  and 
the  unknown  electric  field   E(s,k).   The  coefficients  of  the  moments 
and   E  are  Integrals  which  are  known  In  principle.   The  same 
statements  apply  to   g(s,k,w).   Consequently  all  of  the  moments  of 

/N  ^v      /\      ^ 

g  and   g  are  also  expressed  as  linear  functions  of  n,n,u,u,T,T 


-^  -v^ 


and   E. 

In  the  remainder  of  this  section  we  reduce  the  problem  to  that 
of  solving  for  u,u,  T,T  and   E,   and.  In  so  doing,  construct  the 
dispersion  relation. 

C.   Construction  of  the  Dispersion  Relation 

We  construct  the  moments  In   (s,k)   space  directly  from  the 
definitions  (2.12)  and  the  identity 


r 


g(s,k,w)  = 


\J 


(s,k,neirid^^ 


n(s,k)/n  = 

M^  O 


(s,k,Od^l  =  g(s,k,w) 


(5.10) 


\J 


w=0 
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(s,k)  =  Ug(s,k,Od^^  =  T  ^w  S(s,k,w) 

r 


u 


w=0  ' 


(3.11) 


T(s,k)/T  = 

A-^      O 


vJ 


-^  -    1)  g(s,_k,4)d^e 


'»-s' 


(5c2)"^V^g(^'i^^)  +  g(^'i^^)] 


w=0 


/\  /\ 


The  moments   n,  u,   and  T  are  computed  from  g  In 

A'-i 

exactly  the  same  way.   The  results  for   g  and   g  are  parallel, 
differing  only  In  the  presence  of  the  caret  (  )  over  the 
appropriate  ion  variables.   All  that  we  need  to  complete  the 
system  of  equations  is  a  relation  between   E  and  the  moments  of 
g  and   g.   This  relation  comes  from  the  transformed  Maxwell's 
equations  (5. 12).   The  initial  values  are  not  written  but  implied. 


k-B  = 


0 


k-E  = 


iq/e^  > 


IkxE  =   sB 


'-^    /VV-i 


-ikxB  =   (i,   J   +   sE/c' 


(3.12) 


We  eliminate   B  from  the  equations  (3-12)  and  get  the 
required  relation  between   E  and  the  moments  of   g  and   g. 


-1/,  2_2 


2n-1 


2- 


E  =  e'^Ck'^c^  +  s'^)    (i  k  c  q  -  sj)  . 

The  equation  of  continuity,  which  is  a  consequence  of 
Maxwell's  equations,  is 


sq 


i  k-  J  = 


0 


(3.13) 


(3.14) 
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We  use  equations  {2. 3),    (2.11)  and  (3-1^)   to  write   E  In 


'^^ 


terms  of  u   and   u. 


E  = 


n 


e  s 
o 


s  (u-u)  +  p(p- (u-u) 

2     2 
s   +  p 


(3.15) 


where 


p  =  kc. 


The  vector   p,   which  comes  from  Maxwell's  equations,  is 
treated  differently  from  the  vectors   K  =  kC   and   K  =  kC,   which 
come  from  the  Boltzmann's  equations.   One  reason  for  such  separate 
treatment  is  that  in  some  limits  where   K  and  K  are  small   p 
Itself  could  be  finite  or  large. 

The  evaluation  of  the  moments  given  by  equations  (5-ll) 
simplifies  in  the  rectangular  coordinate  system  in  which  ^        is 
directed  along  one  of  the  axes.   ¥e  represent  any  vector  X  by 
its  three  coordinate  components  as  follows: 


:l^  (x^.x^.Xj) 


In  this  coordinate  system 


B  =  B  (0,0,1) 


With  no  loss  in  generality  we  specify  that   k   lie  in  the 


-"-lo. 


1-3  plane 


k  =  k( sin  6 ,    0,  cos  9  ) 


where  0      is  the  angle  between  the  directions  of  k   and   B 
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The  system  Is  now  completely  symmetric  about  the  B   axis. 
From  this  point  on  we  replace  the  X   and  X   components  of  any 
vector  X  with  the  combination  X   and  X   ,  where 

After  performing  the  operations  indicated  in  equations  (3.II) 
we  get  the  equations  for  the  moments  of   g. 

-r        ,sin  Q    cos  0w    ,     /Sin  9    cos  0w 
""^   ^   ^3^33  ^  ^+  ^ 2— ^V  ^  ^-  ^ 2 ^h- 


(3.16) 


-  iCb  cos  0  I^m  > 

.  2 
u  =  a^(sin  Q    cos  0)l^,  +  a,  I,,  +  a   ( — 7^ — )  I, 

-  iCb  sin  6  I ,  rp   j 

2 
u  =  a  (sin  0  cos  0 )l^_  +  a   ( — ^ — )l   +  a_I__  -  iCb  sin  6    1_^    , 

-,^  T     .       „  T    ,  .    sin  0  T-    ,  .    sin  0  -r     ^-,  -r 
-3C  Tji-  =  la^cos  0  l^r]-,  +  ia_|_  — ^ —  I_|_r]-,  +  ia_  — ^ —  ^_t  "  '"'^-^tT  * 

The  various  integrals   I. .   are  given  in  Appendix  A.   The 
four  equations  (3.I6)  with  the  corresponding  four  equations  for 
the  ion  quantities  and  the  equations  for  E  ,  E   and   E_   given 
by  (3.15)  are  a  closed  set  of  equations. 

When  we  use  equations  (3.I5)  to  eliminate   E  we  end  up  with 

/N        /\        /S 

eight  equations  .for  the  eight  unknowns  u  ,  u_^,  u_,  T,  u.^,  u_^,  u_, 
and  T. 
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The  remaining  substitutions  in  a.,      a,      b  and  b   (given  by 
equations  (2.13))  are 

u  =  (mu  +  mu)/(m+m) 
and  (3.17) 

T  =  (T+T)/2  . 

Equations  (3.I7)  come  from  equations (2. 2 )  and  (2.11 ). 
The  final  form  of  the  equations  is  the  matrix  equation 

AU  =  V  (3.18) 

where   U  is  a  column  vector  with  the  components  in  descending 

,^         /\         y^         /\ 

order:   u^ ,  u  .  u  ,  T/T  ,  u^ .  u  ,  u  ,  T/T   .   The  vector  V 
3   +   -     o   3   +   -     o 

contains  the  initial  values,  and  A   is  an   8x8  matrix  of  the  form 

^2 
A  =  (  \  .  (3.19) 

The  matrices   D  ,  D   D  and  D   are  each   4x4.   Before 
writing  down  these  matrices  we  need  to  define  certain  quantities 
which  appear  in  them. 

M  =  mm/(m+m.)  =  reduced  mass,  (3.20) 

CD  =en/(eM)  =  plasma  frequency  squared, 

K  =  kC  ,    K  =  kC  , 


Vp   =   v^M/m    , 

"2  = 

^   VgM/m   , 

V      =   v-v^      , 

V        = 
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76=  i^) 


2 

CO  s 

o 


m''   2,  2  '   /-    '^    2,  2 
s  +p      ^    I  m  /  s  +p 


2  2  ,22 

n  =  (— J 


,,,      CO  p  ^     A,\   CO  p 


ct  =  V  -  y^       ,   cc  =   V      -    yC 

The  matrices  D   and  D   are  given  below  In  equations  (3,2l) 
and  (3-22)  respectively.   In  them  there  are  the  additional 
Integrals  H  ,  H  ,  H  ,      and  H™  which  are  defined  In  appendix  A 
along  wl.th  the  Integrals  ^-^^>    ^^+'    s^"^-   ^"  order  to  get   D   and 
D   from  D   and  D   respectively  we  need  only  place  carets  above 

/s. 

the  appropriate  quantities.   Thus   D   Is  obtained  from  D   by 
placing  carets  over  a,ri,C,v   and  the  Integrals.   A  caret  over  the 


/\   /\ 


appropriate  Integral,  say   I^^?   means   K,v,n,|-i  and  i/   go  to  K,  v, 


^      /x 


fi,  \i     and  f   In  the  Integral   I_,_,   appearing  in  appendix  A, 
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The  dispersion  relation  Is 

det  A  =  0  .  (3-23) 

The   8x8  determinant  of  A   can  be  reduced  to  a  more  compact   4X>4 
determinant  with  the  use  of  standard  transformations  such  as  may  be 
found  In  A.C.  Altken's  book     .   It  Is,  however,  easier  to  take 
limits  with  det  A  Itself  so  that  our  future  manipulations  start 
from  equation  (^•23). 

The  elements  of  A   are  complicated  functions  of   s,  k,  9, 
the  physical  parameters   0,^2,  m,m,  v-,,v    and  to  ,   the  sound  speeds 
C   and   C,   and  the  speed  of  light   c.   The  variable  k   appears 
In  the  products   K  =  kC,   K  =  kC,   and   p  =  kc.   In  most  of  the 
analysis   p   Is  treated  as  Independent  of   K  and   K.   For  example 
we  may  seek  roots  of  the  dispersion  relation  In  which  K   (and 
hence   K)   Is  prescribed  to  be  real  and  small  with  respect  to 
ffif.   Then  we  find   s   as  a  function  of   K  while   p   Is  an 
Independent  physical  parameter.   On  the  other  hand  If  we  seek  roots 
which  are  characteristic  of  electromagnetic  propagation  we  find   s 
as  a  function  of   p  while   K  and   K  are  treated  as  small 
quantities  with  respect  to  p  by  virtue  of  the  fact  that  the  plasma 
Is  non-relatlvlstlc. 

At  this  point  It  Is  useful  to  sketch  the  connection  between 

[21 
the  lowest  order  results  of  R.  L.  Llboff   -'  and  those  found  by 

using  the  IFF  model.   Llboff  used  a  generalization  of  the  Krook 

model  to  construct  a  set  of  moment  equations  similar  In  kind  to 

equations  (5.16).   He  then  constructed  a  dispersion  relation  and 
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Investigated  two  limits  near  k  =  0,   that  is,   where   k  =  k  times 
a  suitable  characteristic  length  is  small.   In  the  first  limit  he 
sought  roots  (in  our  notation)  in  which   s  =  s^  '^  +  s^  '^k  +  .  .  . 
He  found   s   '^   by  dropping  k  and  higher  order  terms.   In  the 
second  limit  he  sought  those  roots  in  which   s  =  s  ^k  +    . . .    . 

Again,  he  evaluated   s   '^  by  consistently  dropping  higher  order 
terms  in  his  expansions.   To  the  order  of   s  ^      in  the  first  limit 
and   s^  '      in  the  second  limit  the  IPP  model  gives  the  same  results 
as  Liboff 's.   It  is  only  to  higher  order  that  the  two  models  dif- 
fer.  One  may  uncover  all  of  Liboff 's  roots  of  the  first  kind 
(s  =  s^  ^   =   const.)  from  equations  (3-23)  by  suitably  expanding  the 
integrals  and  dropping  terms  of  the  proper  order.   To  get  his  roots 
of  the  second  kind  (s  =  s^  '^k)   one  uses,  instead,  equations  (3-13)j 
(3.16)  and  the  equations  for  n   and   n  which  we  have  not  written 
down  since  we  do  not  need  them  in  the  sequel.   Since  both  theories 
are  equivalent  to  this  order  we  do  not  give  the  details  of  these 
computations. 

The  next  section  contains  the  principal  results  of  this  in- 
vestigation:  the  study  of  the  Larmor  resonances.   It  is  with 
these  resonances  that  the  predictions  of  the  two  theories  become 
different. 
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^.  -The  Larmor  Resonances 

A.   Introductory  Remarks 

The  Larmor  resonances  are  roots  of  the  dispersion  relation 
having  the  form 

s  =  iN|fif  +  R^  (4.1) 

where  N   is  an  integer  and  R^^  is  some  function  of  k,  0,  p,  N 
and  the  physical  constants.   There  is  a  corresponding  set  of 

A. 

resonances  involving  the  ion  Larmor  frequency  n  .   To  be  sure, 
this  is  a  statement  about  dispersion  relations  rather  than  about 
physically  measured  resonances.   For  this  reason  the  nature  of  the 
resonances  depends  on  the  collision  model  from  which  one  derives 
the  dispersion  relation.   Hopefully,  measuring  techniques  will 
improve  to  the  point  where  the  predictions  of  competing  models  will 
be  subject  to  experimental  verification  or  rejection.   Indeed,  the 

predictions  based  on  the  IPP  model  differ  significantly  from  those 

[l  2l 
based  on  the  Krook  model  and  its  modifications  '      . 

We  need  certain  representative  lengths  for  the  rest  of  the 

discussion.   These  lengths,  given  by  equations  (4.l),  are  defined 

up  to  multiplicative  constants  which  do  not  affect  orders  of 

magnitude. 
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A  =  the  wavelength  =  l/k, 

A  =  the  collision  mean  free  path  =  C/v., 

Aq  =  the  electron  Larmor  radius  -   C/\Q\  ,  {^•2.) 

A„  =  the  ion  Larmor  radius  =  C/fi^ 

Ad  -   the  Debye  distance  =  (C  /cd  )'^    . 

We  do  not  distinguish  between  ion  and  electron  mean  free  path 
since  equations  (2.IO)  and  (2.14)   indicate  that  they  are  of  the 
same  order. 

The  condition  that   k  be  real  is  imposed  throughout.   In  the 

region  of  k   space  for  which  ( Ap,/A )  «  1,  ^^c/'^^   «  -'-   ^^^ 
(a  /'a)  «  1   we  find  two  significant  differences  between  the  IFP 
and  the  Krook  model  Larmor  resonances.   The  first  is  with  respect 
to  the  angle  Q      between   k  and  B  .   Prom  equation  (3.23)  one 

'\^^  '^O 

finds  roots  of  the  form  (4.l)  for  arbitrary  N  at  angles  6      in 
the  range 

0  <  0  <  TT  . 

On  the  other  hand  the  Krook  models  predict  roots  for  arbitrary 
N   only  at   0  =  tt/2        and  an  extremely  small  angular  cone  about 
©  =  it/2     which  decreases  as   [n]   increases.   This  point  will  be 
clarified  in  subsection  B. 

The  second  significant  difference  is  in  the  multiplicity  of 
roots  of  the  form  (4.1)  for  each  given  N.   From  the  IPP  dispersion 
relation  we  find  an  infinite  number  of  roots  for  each  N  which 
differ,  roughly,  by  integral  multiples  of   v:   from  the  Krook 
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models  we  find  but  one.   This  reflects  the  fact  that  the  IPP 

r  17 1 

operator  is  unbounded  and  has  a  discrete  spectrum     (zero  and 
the  negative  Integral  multiples  of   v),  while  the  Krook  operator 
Is  bounded  and  has  but  one  eigenvalue  (-v). 

One  of  the  uses  of  a  collision  model  Is  In  the  study  of 
colllslonless  limits.   By  colllslonless  limit  we  mean  a  process 
by  which  the  two  body  collision  frequencies  become  negligible  In 
some  sense.   The  most  Interesting  such  process  Is  non-uniform  in 
that  it  retains  the  concept  of  collision  mean  free  path.   To  reach 

it  we  require  that   (^^/^)  «  1?  ^^c/'^^   "^"^  ''''   ^^^  keep   (A  /a) 
=  K/v   fixed  and  much  smaller  than  unity  as   v   becomes  small.   As 
we  show  in  subsection  C,  for  each  Integer  N  there  will  be  a  large 
number  of  Larmor  resonances  with  negligible  damping  at  arbitrary 
angles  G    in  the  range   0  <  0  <  ir.   The  degeneracy  of  such  resonances 
about  each   s  =  iNln[   increases  as   v   becomes  smaller.   These 
results  reflect  the  properties  which  distinguish  the  IPP  from  the 
Krook  models.   The  latter  models  cannot  give  such  results  for 
the  sdfne  reasons  indicated  above  for  the  case  where  colllsion^^SJ 

A 

are  not  negligible.-      •  ' "  -  - 

In  subsection  D  we  consider  the  uniform  colllslonless  limit 
in  the  range   (Ap^/x)  «  1   and   (A^/a)  «  1.   ¥e  reach  this  limit 
by  letting  the  collision  frequencies  vanish  while  all  other 
quantities  remain  fixed.   In  the  uniform  limit  both  Krook  and  IFF 
models  lead  to  dispersion  relations  which  are  equivalent  to  that 

r  -1  o  "I 

derived  by  I.B.  Bernstein    ^    from  the  colllslonless  Boltzmann 
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equations.   The  only  reason  for  elaborating  on  this  case  is  that 
Bernstein's  results  for  the  Larmor  resonances  were  given  for 
longitudinal  electric  waves.   That  is,  he  applied  the  constraint 
that   E  be  parallel  to  k.   This  considerably  simplifies  the 
dispersion  relation.   R.M.  Lewis  and  J.B.  Keller  -^  ,  in  their 
study  of  a  Krook  model  dispersion  relation,  found  that  this 
constraint  imposes  certain  conditions  which  must  be  satisfied  by 
the  transverse  components  of  the  dispersion  tensor.   The  conditions 
are  not  satisfied,  however,  in  the  range   (^o/^)  «  1  at  6   =tt/2. 
This  is  not  to. say  that  Bernstein's  results  are  wrong  in  this 
range.   Rather,  one  must  resolve  the  question  of  whether  or  not 
the  transverse  electric  waves  are  negligible  by  studying  the 

analytic  properties  of  the  initial  data.   The  dispersion  equation 

* 
itself  does  not  give  answers  to  such  questions  .   For  this  reason 

we  present  the  results  without  the  constraint  to  longitudinal 

waves.   The  principal  effect  of  relaxing  the  constraint  is  that 

there  are  a  threefold  rather  than  a  singlefold  infinity  of  Larmor 

resonances  for  each  species. 

In  subsection  E  we  continue  our  study  of  Larmor  resonances  by 
considering  a  plasma  with  finite  collision  frequencies. 

The  first  study  is  for  9  in  the  range  0  <  0  <  tt  and 
(A^/a)  «  1,  {\/^)  «  1  and  (A  /a)  «  1.  At  0  =  Tr/2  we 
compare  the  IPP  results  with  those  of  R.L.  Liboff    .   ¥e  then 


*  For  this  observation  we  are  indebted  to  Prof.  H.  Weitzner. 
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study  the  fundamental  Larmor  resonances  at  0=0   In  the  limit 
(a  /a)  «  1.   Finally  in  subsection  F  we  study  the  Larmor 
resonances  in  the  range   (A  /a)  »  1  at   6=0. 


B.   Analysis 

Except  for  the  case   (A  /A)  »  1   the  computations  leading  to 
the  Larmor  resonances  rest  on  the  following  analysis. 

We  make  the  change  of  variables   y  =  z/[i7|   in  the  integrals 
defined  in  appendix  A.   It  is  convenient  to  make  all  the  frequencies 
in  the  matrices   D   and  D   (equations  (3-2l)  and  (3.22))  non 
dimensional  with  respect  to   [fi[.   The  dimensionless  frequencies 
carry  a  bar  (  )  :  s,  p,  cd^,  v,  etc.   The  ion  integrals  and 
frequencies  are  treated  identically  with  respect  to  \Q,\       since  we 
focus  our  attention  on  the  electron  resonances.   There  is  a 
parallel  development  leading  to  the  ion  resonances. 

The  elements  of  D   and   D   contain  typical  integrals  of 
the  form 


l{a-^,d^)   =  I  dz  exp|7  -  id^v   +   id^)^^   ,         (4.2) 


o 


where  the  integers   d   and  d   have  the  ranges 


0  <  d-^  <  4    and    0  <  [d^f  <  2  , 

2 


7  =  T  -  (s  +  VM.  )z  +  7  +  ;^_i_  +  ^_  ,  (4.3) 
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T  =  ^L 


3 


2  ^^c  /     2    K^cos^g 
l^j_  cos  if    ,   n      = 


^3 


2  ^  K^sln^e     2 
V  +fi 


2   2 


+ 


-^^  =  -  (v+l(^)/(vtifi)  , 


■^  2  ^-vz 

^+  =  (kL?/2)  exp  )  +  i^  -  (V  ;  i): 


The  quantities  ^^,      ^   and  ^   are  bounded  functions  of   z 
In  the  range  of  integration.   Therefore  the  expansion 


—7  00 

3  - 


(hP 


ri3=0 


^ 


3' 


converges  absolutely  as  a  function  of   z   as  do  the  similar 
expansions  for   e   and   e   .   After  substituting  the  expansions 
into  (4.2)  we  Integrate  term  by  term  for  real  part   (s+vp.  )  >  0. 
The  result  after  a  straightforward  change  of  summation  variables 
is 


I(d^d2)  - 


(4.4) 


+  00 


zszz 


p=0  N=-oo  L=0 


{-^?)^^?/2)l"l-^2L  e"* 


PJ  (  [n|+L).'L.'  U+v(^i^+[N[+2L+P+d^)-i(N-d2) 


For  k   (and  hence  [x^)      real  the  analytic  continuation  of 
l(d  ,d  )   into  the  complex   s   plane  is  direct.   The  function 
l(d  ,d  )   is  analytic  except  at  the  poles  in  the  denominators. 

There  is  but  a  finite  number  of  terms  of  the  triple  sum 
contributing  to  each  pole.   This  number  depends  on  the  value  of 
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2L  +  P.   Since  the  poles  are  distinct,  one  can  show  that  at  an 
arbitrarily  small  distance   e  >  0  from  one  of  them  the  function^ 
I(d-,  ,d  )   satisfies 

That  is,   I(d  ,dp)   is  bounded  away  from  the  poles. 

There  is  an  expansion  similar  to  (4.4)  for  the  ion  Integrals. 
The  ion  poles  are  distinct  from  the  electron  poles  by  virtue  of 
the  mass  differences. 

We  use  the  expansions  (4.4)  to  find  the  Larmor  resonances  as 
roots  of  the  dispersion  equation  (3-23)  for  real   k   in  the 
parameter  range: 

(A^/A)  «  1,    (A^A)  «  1   and  {\/^)   «  1-     (4.5) 


In  terms  of  the  components  of  I-l  and  |j,  the  inequalities  (4.5) 
are  equivalent  to 


H^l  «  1,  \lij    «  1,  Im-^I  «  1  and  \iij    «  1.       (4.6) 


The  poles  of  expansion  (4.4)  suggest  the  possible  Larmor 
resonances.  For  example,  we  fix  N  =  N  >  1  and  guess  that  a 
possible  root  might  be 

i  +  Vm,2  .  i^  +  (l/r)(n2/2)  °(-n;)  °  ,  (4.7) 


40 


where 


s   =  IN   -  vN  , 
o     o      o 


y  and  the  numbers   N   and  N   are  to  be  determined.   We  look 
'  o        o 

for  the  lowest  order  root   in  \x,      a  root  for  which  y  Is 

Independent  of  \i     and  the  magnitude  of  K   and  for  which  the  sum 

I     " 
N  +  N   Is   a  minimum.   Inspection  of  the  Integrals  appearing  In 

(5.23)   shows  that  the  lowest  order  terms  containing  y   come  from 

1(1,-1)  and  the  product  m-x  1(2,-2).  We  substitute  (4.?)  Into  l(l,-l) 

and  find  that  for  N  =  N  -1   and  N  =  0 

00  o 

l(N^-l)^ 

1(1,-1)  =  _  ^   +  y  ^,^   _^s,   +  o(n^)  +  o(-n^)  . 

s  +V-1        ^  o 
o 

Similarly  |j.j_  l(2,-2)   contributes  a  term  In  y  plus  higher 

P        p 
order  terms  In   |i,j_   and   |_l  .   All  of  the  terms  In  (3.23)  of  order 

la-^ ,  \ij_,    M-^,  |_Lj_   and  higher  are  dropped;  thus.  In  the  terms  a,  p, 

Tj   and  In  the  corresponding  Ion  terms   s   Is  replaced  by   s 

In  a  similar  way  we  find  the  other  roots  suggested  by  the 

poles  of  the  Integrals  (4.4).   For  each  such  root  the  expansion  of 

the  determinant  (3-23)  leads  to  an  equation  which  Is  either  linear 

or  quadratic  In  the  unknown  y.   We  then  find  y  and  verify  the 

existence  of  the  Larmor  resonances  as  roots  of  the  dispersion 

relation.   The  same  procedure  with  appropriate  modification  leads 

to  the  Ion  resonances  and  to  the  various  limiting  cases  to  be 

discussed. 
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At  this  point  we  focus  attention  on  the  angular  behavior  of 
the  Larmor  resonances.   We  take  as  typical  the  resonance  described 
by  equation  (4.7).   The  root  may  be  written 

s  +  vn2  =  lN   -NV  +T-\'=|-  ^^}  ,  (4.8) 

where   K^  -  K^/\l^      and   N  >  1. 

The  only  restriction  on  angle  is  that  6      lie  in  the  open 

interval   0  <  0  <Tr  .   The  end  points  present  the  only  non- 

uniformitv  since  for  N  >  1   there  is  no  resonance  at   0=0  and 

o  ' 

Q   =  IT.      Within  the  allowed  angular  range  the  form  of  the  resonance 

changes  continuously  with  angle.   The  Krook  models,  however,  behave 

t 

differently.   To  show  this  we  study  the  Krook  Integral   I    which 

corresponds  to  our  integral   l(l,-l).   The  example  comes  from  the 

fiql 

paper  of  R.M.  Lewis  and  J.B.  Keller  -^  .   In  our  notation 


l'  =  I  dz  exp(^'+iz)  (4.9) 

-^o 

where  _„  ^   _„ 

^     =  -K^  -  (s+v)z  -  -^  +  —  (e^^+e" 


K^  =  K^  cos^e  ,   K^  =  K^  sin^e,   K^  =  K^  +  K^ 


3  '   ^  '        3 


K 


2 

i   +iz 


The  functions  -^  e     are  bounded  in  the  region  of 
integration  so  that  the  expansions  of  the  exponentials  converge. 
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After  some  manipulation 
+00   00 


I  =  e 


2 


2/o\  N  +2L 


N=-oo  L=0 


( K^/2 ) 

(  |N|+L}.'Bi 


00 

r 


\j 


dz  exp 


-(s+v-l(N+l))z--S 


K?  ,2 


(^.10) 


We  seek  the  resonance  characterized  by  N  >  1 .   To  lowest  order 

2 
In  Kj_  ,   the  two  contributing  terms  are   N  =  L  =  0,   and  N  =  N  -1 


and  L  =  0.   The  latter  term  Is  designated   I 


N 


N  -1  r 
,    (K^/2)  °    ' 

\  ~         (N  -1)J 


^ 


dz  exp 
o 


K?  z2 


■(s+v-lN^)z  -  -^2 


(4.11) 


We  require  that   0  <  K?  «  1.   Let 


T  =  I K  I z   and   a 


(s+v-lN  ) 
J 2^ 

IKJ 


(4.12) 


then 


where 


,—0 


(Ky2) 


N  -1 
o 


N 


o   (n^-i)j1k.| 


(a) 


(a) 


r 


dT  exp(-aT  -  T  /2 ) 


\J 


o 


(4.13) 


(4.14) 


One  may  get  an  asymptotic  expansion  of   1(q)   ^o^   I'^I  ^^  1  ^V 
Integrating  by  parts  several  times,  and  a  Taylor  expansion  for 
\o\    «  1   by  expanding   e~  '^      and  Integrating  the  required  number 
of  terms.   These  expansions  In  slightly  different  form  are  found 
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In  reference  [19]- 


For 


a|  »  1 
1 


(a) 


"N^ 


-1-0  V 

o  \    a 


(4.15) 


For   I  a  I  «  1 


(a) 


TT 
2 


a  +  0(a'^) 


(4.16) 


We  suppose  that  \o\    »  1   and  substitute  the  first  term  of 
(4.15)  into  (4.14).   The  result  is 


o   N  -1 
( K^/2 )  ° 

"N  ~   (N  -1).' 
o      o 


We  now  guess  that 


1 


s+v-lN 
o 


N  -1 
■V  +  IN^  +  (1/y  )(k2/2)  ° 


(4.17) 


(4.18) 


To 


lowest  order  in  the  small  quantities  K^      and  K^, 


I  = 


1     ,     Y 
i(N^-l}  +  (N^-1).' 


(4.19) 


The  determination  of  y        from  the  appropriate  Krook 
dispersion  formula  then  follows  immediately.   The  procedure  is 
valid  provided  that   |a|  »  1   or  equivalently  that 


Y 


/  '\N  -1 


1 


»  1. 


K- 
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If  we  suppose  that   |  y  |   is   0(l)   the  Inequality  "becomes 


cot  0|  « 


K, 


2 


K 


2^N  -2 


(4.20) 


Now  suppose  that  \o\    «  1   and  that  we  still  wish  to  find  a 
root  of  the  form  (4.l8).   Substitution  of  (4.l6)  and  (4.l8)  into 


(4.13)  _shows  that  this  is  impossible.   For  example  let 
«  1   then  for 


o|  < 


2 


N  -2 
cot  0|  >  (Kf/2)  ° 


(4.21) 


there  can  be  no  root  of  the  form  (4.l8).   The  region  of  validity 

might  be  expanded  for  particular  values  of  the  physical  constants 

and  N   all  of  which  determine  y   .   None  the  less  the  conclusion 
o 

remains  the  same:   there  will  be  an  angular  region  about   77-/2 
beyond  which  there  can  be  no  resonances  of  the  form  (4.l8).   One 
notes  from  equation  (4.2l)  that  the  angular  region  shrinks  with 
increasing  N  .   No  such  restrictions  appear  in  the  Larmor 
resonances  deduced  from  equation  (3.23)- 

r  -1  o  ] 

I.B.  Bernstein    ^    observed  a  similar  angular  restriction 
on  the  Larmor  resonances  in  the  collisionless  case.   More 
precisely,  setting  7=0  in  equations  (4.9)  to  (4.19)  gives 
roughly  the  same  angular  region  of  existence  for  the  Larmor 
resonances  for  both  the  collisionless  and  the  collisional  Krook 
equations.   Again,  the  region  might  be  expanded  or  contracted 
for  special  values  of  the  constants  which  determine  7  . 
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In  general  there  are  no  roots  of  the  form   s  =  _  1-v    to 
lowest  order.   Rather  one  obtains  from  the  dispersion  relation 
(3.23)  the  same  plasma  electromagnetic  roots  found,  for  example, 
in  R.L.  Llbof f ' s  study    .   His  paper  contains  a  detailed  study  of 
these  roots  and  of  the  special  limiting  cases  leading  to  the  ap- 
propriate Larmor  resonances  which  bear  the  indices   N  =  _  1.   We 
do  not,  therefore,  treat  these  limiting  cases  here.   There  are, 

however,  roots  of  the  form   s  +  vti  =  t  i  -  Lv  +  R  ,   where 

2    2     — ?  —2 
L  =  2,5,4...,  and  R    is  a  function  of  \i^,    \x^,    6,   00  ,  p  ,  L,   the 

frequencies  and  the  masses.   Such  roots  can  not  be  found  from 

dispersion  relations  derived  from  the  Krook  models. 

Certain  comments  are  in  order  before  we  proceed  further.   By 

seeking  the  lowest  order  roots  we  lose  three  other  possible  roots 

at  each  resonance  since  each  pole  appears  linearly  in  each  element 

of  the  matrices  D   and  D   of  equation  (3-21)  and  (3.22).   To 

obtain  the  higher  order  roots  one  would  need  to  keep  enough  orders 

2        2 
of  \i^      and   M-j_   so  that  the  expansions  of  the  elements  of 

equation  (3-23)  would  be  consistent.   Evaluation  of  the  dispersion 

equation  would  be  equivalent  to  solving  a  quartic  equation  in  the 

unknown  function  y  which  now  depends  on  the  components  of  \x    . 

In  order  to  simplify  the  display  of  certain  formulae  involving 

the  electron  resonances  we  consider  the  limit   m  —7  00.   In  this 

limit  the  ions  become  a  uniform  neutralizing  background,  the 

collision  frequencies   v,   and   v„   vanish,  and  the  dispersion 

equation  (3-23)  reduces  to 


he 


det   D  =  0  .  (4.22) 

The  conclusions  drawn  from  (4.22)  relevant  to  the  Larmor 

/\ 

resonances  are  also  valid  for  m  finite.  The  only  difference 

/\ 

appears  in  the  y's   which  for   m  finite  contain  additional  terms 

/\        /\  

involving  multiples  of  the  elements  of   D-,  and   D   and  the  a  and 

rj'  terms  again  contain  the  appropriate  ion  contributions.   The 

elements  of  D   and  D   are  finite  since  the  poles  In  the  ion 

integrals  are  separate  from  the  poles  in  the  electron  Integrals. 

(This  conclusion  is  a  consequence  of  the  analysis  following  equation 

(4.4l).)    Since  we  can  translate  the  results  for  ion  mass  infinite 

we  can  also  translate  statements  about  electron  Larmor  resonances 

into  statements  about  ion  Larmor  resonances. 

C,   The  Non-Uniform  Collisionless  Limit 

We  study  the  non-uniform  limit  by  first  choosing  the  small 
parameter   e  «  1   and  then  specifying 

V  =  0(e) 
and  (4.23) 

K  =  0(6^) 

so  that   (A  /a)  =  (kA)  =  0(e)  . 

The  choice   K  =  0(£  "*"  )   where   A  >  0   would  do  as  well.   It 
is  easier  to  follow  the  orders  of  magnitude  consistently  with 
(4.23).   We  illustrate  the  computations  in  the  limit   m  — >  oo  and 
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assert  that  the  conclusions  are  valid  for  finite   m  by  virtue  of 

the  discussion  at  the  end  of  section  4B. 

A  consequence  of  equations  (4.23)  is  that  \iij_\    is  at  most 

of  order   e   and   ||a^l   is  at  most  of  order   e.   One  may  use  the 

expansions  (4.4)  in  the  dispersion  equation  (4.22).   A  natural 

subdivision  in  this  study  is  given  by  the  angular  domains  defined 

by  (sin  0    cos  0)^0,  cos0  =  0,   and   sinS  =  0. 

(i)  sin0  cose  y^   0 

We  look  first  for  the  root  about 

s"  =  iN   -  VN  ,      N  >  1. 
o      o      0      0 

Our  previous  guess  (equation  4.8)  takes  the  form 

N  -1  i(N  -l)^ 
-   -  2   -     (^^?/2)  °   e    °    N  (i+V-i) 

s  +  v,2  _  3^  = (N^-1):  r '—^ '  (^-2^) 

where  y     is  the  unknown  quantity.   We  substitute  (4.24)  into  the 
integrals  of   D   given  by  (5-2l)  and  drop  all  terms  of  order   £ 
and  higher.   To  lowest  order 

I__  -  H_  =  l/i(N^+l)  ,  (4.25) 

I++  =  H^  =  (l+Y)/i(N^-l)  , 

I^_  =  I^^  -  I^_  =  I^^=I^^  =  I_^  =  I^^  =  0. 


Also   to   lowest   order 


.^2    s    /^^2     ,    ^2 
o      o 

-2-2  1^  f-2     ,    -2 


a  =    -CD      s   /(s      +   p    )    , 
o       00         J"     ''     ' 


Tl     =    O^^p^/Sjs^     +     P      )         . 
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(4.26) 


If  In  addition   (N  e)  «  1   is  Imposed  and  we  drop  terms  of 
order   (N  e)   and  higher,  then   s   may  be  replaced  with   IN 
unless  p  -  N  =  0  for  some  particular  values  of   p  and  N  . 
After  we  substitute  (4.25)  and  (4.26)  with   s"^  =  IN^   Into  (4.22) 
we  need  only  solve  an  equation  linear  In  y,      the  unknown 
quantity.   Therefore  y  Is  known  to  lowest  order  In  terms  of  B, 
N  ,  a  and  'r\.      Should  y  vanish  for  exceptional  values  of  these 
parameters  one  must  compute   y  to  higher  order  In  s.      The 
evaluation  of  y  establishes  that  the  Larmor  resonance 

"i  =  IN   -  N  V  +  0(£^)  (4.27) 

0      0 

Is  a  root  of  the  dispersion  equation  (4.22).  In  a  similar  way  one 
can  demonstrate  by  the  same  procedure  that  there  are  either  one  or 
two  lowest  order  roots  for  each  L  >  1   of  the  form 

"i  =  IN^  -  7(|N^|  +  L)  +  0(6^)  (4.28) 

where 

In  I  >  1 
o 

and 

L  =  1,2, ... 

The  number  Is  two  If   0([i^)  =  0(e'  )  =  0(m.j_)   and  one  If 

O(n^)  =  0(6^),   O(ii^)  =  O(e^)  .   The  terms   0(6^)   In  both  (4.2?) 

_  2 
and  (4.28)  are  contributed  by  v\x    .   The  terms  Involving  y  are 

generally  of  higher  order  although,  as  mentioned  above,  peculiar 

combinations  of   p^,  co^,  N^,  and  0   might  make  y  =  0(e  ).   We 
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assume  that  In  general  such  Is  not  the  case.   One  may  go  to  the 

limit  £  =  0   in  the  following  way.   For  the  root  (4.27)^  fo^^ 

example,  we  substitute  equation  (4,24)  into  the  integrals  of  (4.22) 

and  then  take  the  limit   e  =  0  and  hence   v  =  0  and  K  =  0  by 

(4,23).   The  elements  of  dispersion  equation  for  y     remain  finite 

and  we  say  that  in  the  collisionless  limit  the  dispersion  equation 

is  compatible  with  the  Larmor  resonance   s  =  IN  .   If  we  fix  N 
■^  0  0 

and  substitute  into  (4„22)  successively  for  each  value  of  L  in 
(4,28)  the  form  corresponding  to  (4.24)  and  go  to  the  limit   e  =  0, 
then   s  =  IN   is  again  a  resonance  compatible  with  the  dispersion 
formula  (4.22).   In  this  sense  the  Larmor  resonance   s  =  IN  ,  for 
N  [  >  1,   is  infinitely  degenerate.   For  N  =  t  1   the  limit 

O  '  JO  Q 

e  =  0  does  not  exist  except  when  cd  /"v   remains  finite  as   v 

0/ 

goes  to  zero.  We  see  this  by  a  closer  Inspection  of  the  term  I  . 
For  example  suppose  we  look  at  the  resonance  about  s  +  v|-l  ,  =  i-2v. 
To  lowest  order 

I^^  =  ('s+V^L^-i+V)~^  +  {-\i^){'s+v\i^-±+2v)~-^      .  (4.29) 


If  we  suppose 


s  +  vu   -  i  +  2v  =  (-m-^)/y  (^.30) 

then  to  lowest  order 

I++  =  -  i  +  7  .  (^.31) 

V 

From  (4.31)  Y  =   v~   and  as  a  next  approximation  one  might  write 
y  =  (1+vA)/~  ,  substitute  into  the  dispersion  relation  (4.22),  and 
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find  that  H   and  the  term  I    give  with  I    a  quadratic 
expression  for  A.   If  go   remains  finite  as   e   and  hence  "v 
vanish,  the  terms  In   I    blow  up.   However  If  <^p, /v   remains 


++ 


o 


finite  as   v   goes  to  zero,  and  If  we  postulate  that   (od  y) 
remains  finite,  we  conclude  that  the  Larmor  resonance   s  =  1   Is 
consistent  with  the  dispersion  equation. 

In  general  for  0  <  e  «  1   one  can  conclude  that  there  are 
resonances  of  the  form 


t  1  -  v(l+L)  +  0(£^)  , 


(^.32) 


where   L  =  1,2,3  • • • 

but  that  the  limit   e  =  0  does  not  exist  in  the  sense  described 

above  for  |N  f  >  1. 
o 

(11,  )  COS  e  =  0 

In  this  case   fj,^  =  0  and   [j,j_  =  [x  in  the  Integrals  of 

—2  — ? 
Appendix  A.   Assuming  that  no  combinations  of  oo  ,  p   and   N   will 

cause  the  lowest  order  determination  of  the   y  s  to  vanish  we  may 

write  the  dispersion  equation  (4.22)  in  the  form 


where 


det  D^  = 


(-1+a  I^^)  det  D-, 

DD  1 


0 


(-l+(a  -  Ti/2)H^) 


+  ■ 


(-  Ti/2  H  ) 


(-  Ti/2  H^) 


(_l+(ct  -  ii/2)H  ) 


(^.33) 


(4.34) 


In  the  exceptional  cases  one  might  have  to  restore   a  I 


+- 
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to  (4.54)  and  include  the  row  and  column  containing   l_^rp  ,  I_rp  , 
Hrp  and   Imm  .   We  suppose  that  the  exceptional  case^do  not  arise. 

Of  the  lowest  order  roots  (4.28)  those  for  L   odd  no  longer 
appear  at   cos9  =  0  when  we  follow  the  usual  procedure  for  finding 
the  roots  of  (4.34).   However,  at  angles  arbitrarily  close  to  ti/2 
equation  (4.22)  does  give  roots  for  odd  L.   Such  apparent  non- 
uniformities  are  a  result  of  working  only  with  the  dispersion 
equation  and  ignoring  the  initial  data.   A  more  obvious  example 
occurs  at   sin0  =  0. 

The  term   (-1+a  I^^ )   in  equation  (4.33)  does  contribute  roots 

for  odd   L.   When   |a^  =  0(8  )   these  were  higher  order  roots  and  did 

not  appear  in  the  cos0  •  sin0  ^   0  case.   When   0(i_l^)  =  0(ii,_j_) 

=  0(e)   these  were  half  of  the  two  sets  of  roots  found  for  each 

L  >  1  .    Equations  (4.2?)  and  (4.28)  remain  the  same  except  that 

0(e  )   becomes   0(£^)   with  the  exception  of  the  combination 

In  I  =  2   and  L  =  0  for  which  case  the  lowest  order  in   e   is 
o 

/    4  N 

0(e  ).   The  limit   e  =  0  does  not  change.   Also  the  resonances 
(4.32)  remain  the  same  as  does  the  conclusion  that  the  limit   e  =  0 
does  not  exist  in  general. 

(iii)   sin0  ^  0 

Here   p.j_  =  0   and   |i.^  =  [x   in  the  integrals  of  appendix  A. 
The  only  terms  contributing  a  Larmor  resonance  are 


(-1  +  a  I^^)  (-1  +  a  I__)  =  0.  (4.35) 
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For  N  =  +1   the  integral   I^  ,  =  \  e^"^'^-""" '^dz  contributes, 


The  roots  are  given  by 

2 

o    L 


B   +   v^^   =    s^   +    (ff-)  y-^  (4.36) 


where 


L-  l,2j3j  •••  J 

Sq  =  :^  1  -  v(L+i)  , 

—2  —2 

,      1     (^°"P  ) 

r     _   _2  — 

Lv   oj   s 
o   o 

We  see  again  that   e  =  0  Is  an  Improper  limit  unless  cd  /v 

remains  finite  as  v  =  0(e)   goes  to  zero  and  co  v  also  remains 

—2  /-  — 

finite.   The  condition  that  to  /v   be  finite  as   v   goes  to  zero 

o  ^ 

can  be  achieved  by  letting  the  electron  density  n    go  to  zero 
since   both  cd   and   v   are  proportional  to  n  .    The  are  no 
resonances  for   [N  f  >  2   at   sin  0=0.   At  angles  arbitrarily 
close  to   0  or  TT  we  get  all  of  the  resonances  (4.28).   As 
before,  such  non  uniform  behavior  is  due  to  the  fact  that  we 
omit  the  initial  conditions  and  deal  with  the  dispersion  equation 
alone. 


D.   The  Uniform  Colllsionless  Limit  at   0  =  ir/2 

Again  we  Illustrate  with  the  dispersion  relation  (4.22)  for 
ion  mass  infinite.   To  reach  the  uniform  colllsionless  limit  we 
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let   V   go  to  zero  and  keep  all  other  quantities  fixed.   At 
9   -  tt/2      the  dispersion  equation  (^.22)  becomes 


det  D   (O)   = 


[1+a   I^^] 


0 


0 


[-1+a  I^^-  I  H^] 


0 


0  i[aI+_-   nH_] 


i[aI^_-TlH^] 


[-l+al_    -  n  H_] 
2 


where 


33 


<j 


,^ 


dz. 


'++ 


00 

-   [e^  fe+i^  +  (K2/2)(l-e+i^)4   dz, 

00 

:       e?    L-1^+(kV2){1- 

o 


-IZ    ^2    1  -, 

)    ?   dz. 


00 

^_=   -K^    re^(l-e+^^)(l-e-^^)dz. 


o 

r 


H      =   + 

+ 


IS 


^^(l-e-^^)dz. 


^  =    _k2_    sz   +   k2/2    (e+^^   +   e"^^)    . 

Typical  Integrals  have  the  form 

-r/ J  \    r  ^  +ldz 
1(d)  =  Je^ 


=  0, 


(^•37) 


(^.38) 


(4.39) 
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for  d   an  integer  in  the  range 


2  <  d  <  2. 


We  may  expand  the  exponential  of  the  bounded  functions 


'T?2 


(KV2)e 


±lz 


and  Integrate  for  real  part   s  >  0.   The  result   Is 


+00 


1(d) 


-k2 


(K^/2)'^'+^^ 


N=-oo  L=0  (  lN|+L)iL' (s-i(N+d)) 


(4.40) 


To  effect  the  analytic  continuation  we  note  that   l(d)   is 
bounded  away  from  the  well  separated  infinite  set  of  poles.   If   s 
is  a  distance   0  <  5  <  |-  from  one  of  them  then 


[l(d)I  <  1/5  . 

Consequently,  since  the  ion  poles  and  the  electron  poles  are 
distinct,  one  may  apply  the  results  for  m  infinite  to  m  finite 
as  before. 

We  suppose  that   K(  =  K/|fi| )  «  1   and  apply  the  procedures 
of  4B  to  (4.37)  and  (4.38).   The  integrals  I_^^  =  H^  ,  I__  =  H_ 
contribute  the  lowest  order  resonances. 


.  =  lN/l+^(K2/2)i"l-li|„| 


(h.m) 


where 


N|  >  1   and 


Ri 


N 


(|N|^2/(^^|jjf))^  (-2/2) 


N^ 


(l-  Z^2/(N2_l)^)  _  N^A-  Z^2/(n2:|+  ^-J/(n2_i)^ 


(4.42) 


-  55  - 


There  are  In  principle  three  roots  obtained  "by  Including 
higher  order  terms  In  K  .  According  to  (4.4l)  there  are  an 
Infinite  niomber  of  Larmor  resonances  at   6  =  7r/2. 

Roughly  speaking,  as  soon  as 

the  resonance  (4.4l)  Is  no  longer  a  solution  of  the  dispersion 
equation.   This  follows  from  the  statements  subsequent  to  equation 
(4.21)  with  the  refinement  afforded  by  the  substitution  of 

(1/r  )  =^4^    . 

We  observe  that  the  formal  limit   |pl  — ^  oo  in  (4.42)  gives 
the  results  one  can  compute  from  I.B.  Bernstein's  longitudinal 
electric  wave  dispersion  relation 


E.   Larmor  Resonances  with  Collisions 

We  choose  the  small  parameter  e  «  1   and  require  that 

K  =  (A^A)  =  0(e),  (4.44) 

v  =  0(1). 
This  choice  Insures  that 


(A^/A)  =  0(e), 

[(1^1   Is  0(e)   or  higher  (4.45) 

f(Xj_|   Is   0(e)   or  higher  . 
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The  order  of   ||a^|   depends  on   cos  Q      and  that  of   |  |j.j_  |   on 
sin  9 . 

Following  the  procedures  of  section  4B  with  respect  to  the 
dispersion  formula  (4.22)  for  m  infinite  we  compute  the  Larmor 
resonances  in  the  ranges:  (sin  6    cos  Q)   ^  0,      cos  0=0,   and 
sin  0=0. 

(i)  (sin  0  cos  0)7-^0 

The  formulas  in  this  range  are  given  for  both   |  [x^  |   and   |  |a,j_  I 
of  order  e  .   There  are  in  principle  four  roots  for  each  resonance. 
As  one  approaches  cos  0=0   or   sin  0=0  the  lowest  order  root 
might  switch  from  one  of  the  four  to  another.   The  order  of 
magnin\»fede  of  the  term  multiplying  the  y's  might  change  but  the 
existence  of  the  roots  will  not.   In  the  mid  range  |  ^.^  [  =  0(e)  and 
[  |i.j_  I  =  0(e)   we  get  the  Larmor  resonances 

s  =  IN  -  v(n2+|Nl+L)  +  0(k2([n[+L-1)^  (^^^g) 


for 


N  /  0,  L  =  O  or  1,  and  ( 1n|+L)  >  1  ; 

s  =  iN  -  V(n2+|N[+L)  +  0(K2(tN[+L-2))  (^^^^) 


for 


N  7^  0  ,  L  =  2,3,  ..  •   . 

There  are  roots  for  N  =  0  and  L  =  2,3,  •..,  but   they 
are  not  Larmor  resonances.   The  resonances  (4.47)  are  contributed 
by  the  integral   Imm  and  represent  therefore  the  influence  of  the 
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the  temperature  terms.   The  resonances  (4.46)  come  from  those 
terms  involving  E  and  the  macroscopic  velocity  u.   If  one  could 
neglect  temperature  by  arranging  the  Initial  conditions  In  some 
suitable  manner,  then  the  resonances  (4.47)  would  come  from  the 
terms  Involving  E  and  u  and  the  exponents  of  the  order  of  K 
would  be  2(|n1+L-1)   as  In  (4.46). 

(11)   cos  0  =  0.. 

Here   p.^  =  0  and  |af  =  KV(  1+v  ).  At  this  point  one  can 

[2l 
compare  the  Larmor  resonances  computed  from  R.L.  Llboff's    Krook 

model  with  those  computed  from  equation  (4.22).   Prom  (4.22)  one 

gets 

7  =  IN  -  V(n?+|N|)  +  0(^2(1^'-!))  (4.48) 

for   In[  >  1.   Prom  Llboff's  model  one  gets 

I  .=  IN  -  V  +  o(k2(|n|-i))  (^^^5) 

for   |n[  >  1.   As   fN|   gets  larger  the  damping  of  the  IPP 
resonances  Increases  while  the  Krook  resonances  have  the  same 
damping  factor.   In  addition  to  (4.48)  we  have  for  each  N  the 
Infinite  set  of  resonances  from  (4.22) 

s  =  IN  -  v(m-x+|nI+L)  +  Oiuif^)  (4.50) 

where 
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A^  =  |N|  +  L  -  1     for  L  =  odd  Integer, 

A^  =  |n|  +  L  -  2     for  L  =  even  Integer,  >  0  . 

Equation  (^.50)  also  gives  Larmor  resonances  for  N  =  1  1  . 
These  do  not  appear  In  the  Krook  theories. 

One  observes  that  In  the  uniform  limit   v   goes  to  zero  both 
equations  (4.48)  and  (4.49)  lead  to  the  resonances  (4.4l).   The 
resonances  (4.50)  do  not  exist  In  the  uniform  colllslonless  limit 
so  that  no  degeneracy  appears. 

The  Infinite  set  of  Krook  resonances  (4.49)  exist  only  at 
cos  0  =  0.   As  one  departs  from  0  =  ±  77"/2   the  resonances  begin 
to  disappear  roughly  at  the  angle  given  by  (4.2l).   This  of 
course  Is  markedly  different  from  the  behavior  of  the  IPP 
resonances. 

(ill)  sin  0=0 

Here   ^ij_  =  0  and   [  p.  f  =  (KA)  =  (^  /^)  =  0(e)  . 

J?  c 

The  Larmor  resonances  come  from  the  two  uncoupled  factors  of 
(4.22) 

-  1  +  al^^  =  0, 

(^.51) 

-  1  +  al  _  =  0  . 

The  results  are 

s  =  t  1  -  v(n2+i+L)  +  0{^f)  (4.52) 
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for 

These  roots  do  not  appear  In  the  Krook  model  theories. 

P.   Larmor  resonances  at   0=0  for   (A  /A)  »  1. 

Here  \x^  =   0     and   ^L  =  ^l^  =  (K/V)  =  {^^/^)   »  1-   At   0  =  0 
equations  (4,.  51)  involve  the  Larmor  resonances.   The  first  of 
equations  (4.51)  is 


CD  s 

1  =  fv  -    ° 


n^^  \  r 


-2  —2 
s'^+p'^ 


2\    .  \„    .       2/-,  „-vz 


exp[-(s+v(l+M,'^)-i)z  +  ^x'^(l-e  ^  )]dz  .       (4.53) 
o 


We  prescribe   |j.  to  be  real  and  positive  and  require  that   v 
be  non  vanishing  although  it  may  be  small.   The  integral  in  equation 
(4.53)  is  related  to  the  incomplete  gamma  function  .   To  show  this 
one  changes  variables  according  to  the  scheme 

t  =  ^2e-vy 


a  =  (s-i+v(l+^i^;)/ V  (4.54) 

a  =  o/v  =  (1  -  cD^Vv  (F^+p^ ) ) 


Equation  (4.53)  transforms  to 

2 

1  =  a  le^^V^'^ft.t^^-l^e-^ 

v^  o 


*  We  are  indebted  to  Prof.  H.  Weitzner  for  this  observation. 
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or 


where 


1   =   a    ]e^   li-^'Tio)    -    Ho,ix^) 


Pio) 


00 

I-'-' 


^  e~^  dt, 


(4.55) 


(4.56) 


J,/^   2>,     u^   -2a I    .(a-1 

J  2 


^^-1)  e-t  dt, 


(4.57) 


The  Integrals  (4.56)  and  (4.57)  converge  for  real  part  of 
a  >  0.   ¥e  will  look  for  roots  of  (4.55)  for  real  part  of   cr  <  0, 
The  analytical  continuation  of  [^{o)      from  real  part   a  >  0  to 
real  part   a  <  0   Is  effected  by  the  Identity 


Fio)  = 


IT 


P(l-a)sin  iro 


(4.58) 


It  Is  convenient  to  Introduce  the  change  of  variables 


Z  =  1-0 


(4.59) 


and  to  consider  real  part   Z  >  0  as  the  domain  of  interest.   With 

(4.59)  and  the  variable  change   t  =  (t/|j,   -  l)   the  integral 

(4.57)  transforms  to 

00 

r 


5(z,n^) 


e-^V^(l+t')-^  dt'  . 
o 


(4.60) 


All  of  the  formulas  and  expansions  relating  to  the  Gamma 
function  are  taken  from  Chapter  XII  of  reference  [20]. 
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Not  only  is   ^  convergent  for  real  part   Z  >  0  but  one  can 
also  bound  It  from  above  and  below.   Prom  above  one  sets   t  =  0 
in  the  exponential  and  gets 

I  ^  I  <  J-   ^  where   Z  =  real  part   Z  >  0  ,         (^.6l) 

or  one  sets   t  =  0   in   (1+t  )      and  gets 

I  M  <  \   .  (^.62) 

To  bound  from  below  one  uses  the  inequality  ( 1+t  )  <  e    and 
gets  the  bound 

-^     <  [  H  .  (^.63) 

Z  +\if^ 

The  dispersion  equation  (^.55)  is  now 

.M-^  2(Z-1) 

^   =  a{^-^ ^-^V   .  (4.6^) 

P  (z  )sin  ttZ 

The  roots  of  (^.64)  are  found  by  a  process  of  successive 
approximations.   (i)  We  need  the  first  approximation  of   Z   to  be 

2  n(  V _-\  ] 

large  enough  to  cancel  the  large  quantity   (e^  |j.  ^    ^).       (ii) 
With  this  first  approximation,  labeled  Z-,,  we  use  (4.59)  and 
(4.54)  to  evaluate  a  ,  and  estimate  the  bounds  of   ^  .   (iii)   The 
next  approximation  is  gotten  by  looking  for  the  root   Z  =  (Z  +h), 
where  h  is  bounded  and  independent  of  [i,      to  satisfy  the  rest 
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of  the  dispersion  formula  (4.64). 

(l)   Since  we  seek  those  roots  for  which  the  real  part  of   Z 
Is  large  and  positive  one  Is  justified  In  replacing   P(Z)  In 
(4.64)  with  Its  asymptotic  expansion 

r(Z)'w  (27r)2  Z^^-i)  e"^  (l+0(|)  +  ..  )  .  (4.65) 

Neglecting  the  higher  order  terms  we  get 

1  =  a  ^  exp[n^+Z-log  ^  -  log  (Z/^^)-(Z-  ^)]  _  J  I       ^^^^^^ 
sin  ttZ-  {2/v)^  { 

We  guess  that 

Z^  =  d^^x^  +  a^log  n   .  (4.67) 

Substituting  (4.67)  Into  the  highest  order  terms  In  the 
exponential  gives 

(l+d   -d   log  d    )iJ.     -    (l+a  log  d    )    log  ij.  =   0    . 
^        000  00 

Thus   d   and   a   are  Independent  of   p.  for   p.  large  and  are 


given  by 


d^  =  5.591, 


a^  =  -  1/log  d^  =  -.782, 


(4.68) 


Therefore   Z   Is  known  precisely. 

(11)   The  first  approximation  to   s   Is 
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s^=  1  -  v(l+d^)^i^,  (^.69) 

where  the  log  [i     term  has  been  dropped  since  it  is  dominated  by 

2 
the  [I     terms. 

The  first  approximation  to  a  Is  denoted  by   cu^  .      We  suppose 
the  |i,^   Is  sufficiently  large  so  that  I «(]_)!   ^s   0(l)   and  write 

^(1)  ^  l^(i)|e^Po  .  (4.70) 

(ill)   The  term  ^  Is  bounded  In  absolute  magnitude  by  l/p.   so 
that  It  may  be  dropped  when  \cl-A      Is   0(l). 
Since   Z   is  known  we  write 

Z^  =  2N^  +  A  (^.71) 


where 


N   is  a  positive  Integer  and   0  <  A  <  2 
The  next  approximation  is 


z  =  z^  +  h  (^.72) 


where  h  Is  to  be  small  compared  with  log  \i.      We  substitute 


{k.JO),    (^.71)  and  (^.72)  into  (^.66)  and  write  _^-|  |a  |   as   e 


and 

h  =  X  +  iy   .  (^.73) 

Then    (^.66)  becomes 

Sin  7r(x+A+ly)   =   e^^'^^   "^^^^S  ^o)   ^^iPo-J  log  do)    . 


We   equate   real   parts   and   Imaginary  parts. 

sin  7r(x+A)e        "^    ~'^'      &     o-i    =   sln(ylogd   +  p    )    cosh  Try, 

Ffi    r      ill        H    1  •  (^-7^) 

cos  7r(x+A)e"^      ^'^"s^^^S  aoJ   ^   cos(ylogd   + -p    )    sinh  Try, 


where 


=    (7r/2)    - 


p'    =    U/2j    -    p 


o 


The  system  (4.74)  actually  has  an  infinite  number  of  solutions, 
To  find  a  few  explicitly  we  look  for  a  solution  for  which  Try  >  5 
so  that  tanh  Try  =  1.   Thus  division  of  the  first  of  (4.74)  by  the 
second  gives  us 

y  logd   +  p   =  Tr(x+A)  +  2Trr  (4.75) 

where   r   is  an   Integer.   Then  we  substitute  cosh  Try  =  sinh  y 
=  e^/2      for  r   large  enough  and  get 

Negative  y's  are  also  solutions.   We  need  only  one  of  the 
solutions  to  exhibit  the  root  to  the  dispersion  equation.   Finally 

i  =  +1  -V(  (1+d  )n^  -  log  n/log  d^  +  (x  +ly)),..  -       (4.77) 
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Prom  the  second  of  equations  (4.5l)  one  gets 

s  =  -i  -v((l+d^)M.^  -  logiVlog  d^  +  (x  +(y  ))  .         (^.78) 

We  cannot  say  that  there  are  an  infinite  number  of  these  resonances 

although  there  are  an  Infinite  number  of  solutions  to  (4.74).   The 

reason  Is  that   |h|   must  be  smaller  than  log  |  |j,[  for  the 
analysis  to  work. 

,/vr  I  2 

Should  I  ct-j  I   be  of  order   |j,   one  may  no  longer  omit  the  term 
In  {)      since  we  know  from  equations  (4o6l)  to  (4.63)  that  to 
lowest  order 

[  H  -  -^    ,  (^.79) 

where   for      [[if    »   1 

(1   +   d)~^  <    0<   d 


o."<©<^^'' 


The  analysis  given  here  Is  good  for   Ola,-,J  <  0(|j.  )   so  that 
for   ||i.[   sufficiently  large  one  may  Ignore  the  contribution  of 

If  we  keep  {[if   fixed  and  large  as   v  — >  0  there  is  one  ; 

limiting  root  consistent  with  the  dispersion  equation  (4.53). 

_2  — 

This  root  appears  when  cu  — ^  0  with   v   in  such  a  way  that 

2/-  °         - 

CO  /v   remains  fixed.   The  root  is   s  =  1,   and  correspondingly 

there  is  a  root   s  =  -1   from  the  second  of  equations  (4.51).   In 
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this  limit   K  also  goes  to  zero  so  that  In  effect  we  have 
computed  one  of  the  zero  order  resonances  discussed  by  Llboff 
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5.    Electromagnetic  Propagation  and  Plasma  Oscillations 

A.   Introductory  Comments 

When  k   Is  parallel  to  the  applied  magnetic  field  B   the 
dispersion  equation  (3.23)  simplifies  after  a  simple  rearrangement 
of  rows  and  col/umns  to 


D  -D  °D^  =  0, 
+   -   3 


(5.1) 


In  (5-1) 


D+ 


(-1+  al_^^) 


++ 


PI 


+  + 


(-1+aI   ) 
++ 


(5.2) 


and 


"3  = 


D. 


31 


D 


32 


D 


32 


D 


31 


(5.3) 


where 


D 


31 


(-l+(a-Ti)H^) 


T 


(-lC(v-iv2)H^) 


(l(a-ri)H^)        C(3-(v-iv„)l 


2^  TT 


and 


D 


32 


(P+Tl)H 


3 


l(p+ri)H^ 


1  2  ^2  H^ 


C     ^ 

?   ?  TT 
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The  determinants   D^^   and  D^^   are  obtained  from  D-,.,   and 
B  by  placing  carets  over  the  appropriate  quantities. 

In  section  5B  we  set  D  =  0  and  compute  the  effect  of  a 

colllslonal  plasma  In  a  magnetic  field  on  the  vacuum  electromagnetic 

?   2 
roots   s  +p  =0.   To  do  this  we  suppose  that   p  =  kc   Is  larger 

than  any  frequency  appearing  In  D   .   We  then  suppose  that 

s  =  lp(l+A)   where   A  Is    complex  and   |a|  «  1.   All  of  this 

enables  one  to  expand  the  Integrals  and  the  quantities   a,  a,    t],    r\ 

In  terms  of  the  frequencies  made  non-dlmenslonal  with  respect  to 

[p[.   The  Ion  mass  remains  finite  In  this  computation. 

In  section  SC  We  set  D^  =  0  and  compute  the  colllslonal 

correction  to  the  plasma  oscillation  mode.   Here  the  plasma 

frequency  cu   is  chosen  larger  than  any  other  frequency  appearing 

In  D   and  again  we  suppose  that   s  =  Ico  (l+A)   where  A  is 

complex  and   Ia|  «  1.   In  the  parameter  range  we  have  chosen  the 

result  for  finite  ion  mass  does  not  reduce  when  we  set  m  — ^  00  to 

[61 
the  result  of  Lenard  and  Bernstein'-   .   We  give  the  reason  for  this 

disparity  at  the  appropriate  point.   At  the  same  time  we  compare 

results  with  those  of  the  Krook  model  theories  of  iLewis  and  Keller 

and  of  Liboff. 


B.   Electromagnetic  Propagation 

The  principal  reason  for  performing  this  calculation  was  to 
check  the  applicability  of  a  formula  describing  the  Faraday 
rotation  of  a  plane  wave  by  a  plasma  in  a  steady  magnetic  field. 
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The  formula,  equation  (5.6)  below,  was  derived  for  propagation 
parallel  to  the  steady  magnetic  field  In  a  colllslonless  plasma 
from  simplified  macroscopic  equations  of  motion.   The  experimental 
arrangement  of  Dr.  Fred  L.  Rlbe  ,    however,  corresponds  to  the 
following  parametrlzatlon.   All  frequencies  are  made  non-dlmenslonal 
with  respect  to  p  =  kc,   and  since  his  source  Is  an  optical  maser 
we  take  a  representative  wavelength  to  be  of  the  order  of  6X10  -^  cm. 
The  small  parameter  for  this  experimental  set  up  Is 

6  -  (m/m)  «  1. 
The  collision  frequencies  as  computed  from  equations  (2.14)  are 


e,  =  -^  =  0(6) 
1   p     ^  ' 


In  addition 


^1  ^ 

l^  -   0(55/2), 

^2   = 

v/p  =    0(5^), 

/\                A. 

e,    = 

e^+e^      ,      e   =   ^l'^^2 

)n 

T 

0 

co2 

^=   0(6)    , 

2p^ 

k2 

%=   0(6)      , 
c 

(5.4) 


(5.5) 


*  Private  communication. 
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and 


0(6) 


A  plane  wave  traveling  in  the  direction  of  B   may  be 
resolved  into  two  circularly  polarized  waves  both  of  which  satisfy 
the  relation   s  =  Ip  in  the  absence  of  plasma.   However,  according 

ri4'i 

to  the  macroscopic  analysis  of,  for  example,  L.  Spitzer     in  the 
presence  of  plasma  the  two  circularly  polarized  waves  denoted  by 


s   and   s    satisfy 


T 


o  +  lil 


.  ^=lP[l«o  (1  -^i^  )J  •  (5.6) 

To  get  (5.6)  from  Spitzer 's  equation  (4-1^)  one  expands  his  formula 
and  keeps  terms  at  most   0(6  ).   The  difference  in  frequency 
between  the  two  polarized  waves 


^^^=.P.„M  (5.,) 

is  a  possible  diagnostic  tool  for  measuring  plasma  densities  since 
p  and   Inl   are  known  and  go   is  proportional  to  n  •   The 
question  to  be  answered  is:   Do  the  formulas  given  in  equation  (5«6) 
and  (5.7)  change  for  the  experimental  parameters  used  by  Dr.  Ribe. 
These  parameters  are  given  in  rough  form  by  equations  (5- ■4)  and 
(5.5).   Prom  D  =  0  we  find   s   and  from  D_  =  0  we  find   s_  . 
For   s   the  dispersion  equation  can  be  written 
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D^=    0  = 


(-1   -    ^P\.^ 


(-1  +  ^^K. 


(5.8) 


where 


) , 


(f )   =   e    -    (|)    , 


P  m     P 


(^)    , 

P 

p 

.      ^\ 

^0            £     M 
0 

> 

/\ 

M  -  -^ 

• 

m+m 


++ 


dt   exp[' 
lo 


^P  P 


SIX  t   +   IX   (1-e         )J 


(5.9/) 


where 


2       K"". 
|i     =  -p     •  - 

P  e 


We  obtain   I    from  I,  ,   by  replacing   (e^  Q,    (i  )  by 


++ 


++ 


/\  /s     /\^ 


{e,   Q.,    \i    ).      The  integral  (5.9)   and  the  dispersion  relation  (5«8) 
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are  the  result  of  changing  the  variable  of  Integration  according 
to  y  =  t/p. 

To  find  the  roots  we  guess  that 


I-  =  [1  +  T  (l+h  +  Ihy)] 

ip  O^    X      ^ 


(5-10; 


where  h   and  h   are  the  unknown  functions  to  be  determined. 
X        y 

We  suppose  that  h   and.  h   are   0(5)   and  verify  the  guess  at 

X       y 

the  end  of  the  computation.   One  substitutes  (5.IO)  Into  the 
formula  for   (— )   given  after  equation  (5.8).   Then  by  the  re- 


maining equations   (  — ),  (H.)   and   (-)   are  known.   The  result  for 
(— )   after  expansion  in  the  small  parameters  is  to  first  order 


(I)  =   5  [(£,-h  )  +  K-l  +  h   -  :5^  )]  . 


■p'    m   '  2   y 


(5.11) 


The  integral   I,,   given  by  (5-9)  becomes 


++ 


I    ^  I  g-l[l+To(l+h^)+  n/p  -l(e-Tohy)]t  exp[-£n2t+n2(l-e-^^)] 
++ 


(5.12) 


After  first  expanding  the  argument  of  the  second  exponential 
to  order   5^   and  then  expanding  the  exponential  Itself  one  gets 

-l(l+To+hx1'+  ^h   -   i(e-Tohy|t^ 


I^_j_  =  ^dt  e 
'o 


20^      6c^ 


(5.15) 
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We  then  integrate   I    according  to  the  formula 

e"^^  t"  dt  =  ^^     for  real  part   a  >  0 

vj  a 

o 

which  is  satisfied  if   e  -  t  h  >  0  as  was  guessed.   The  next  step 
is  to  expand  the  terms 

TT \        for  N  =  1   and  3   to 

lowest  order  and  collect  terms.   We  substitute  the  results  for   !,,, 

I  )    {~) }    (  )>  (~)   a-^d   (— )   into  (5.8),  expand  the  determinant, 
++   P    P    P         P 

and  equate  the  real  and  the  imaginary  parts  separately  and  obtain 
to  lowest  order 

X      2  +  ^2    p  l^-J-^J 

and 

.2 


n     v-2  ^  ^   2   1^ 


Prom  equation  D  =0  we  find  that   h   is  the  same  but  that 

hx  =  -^^4-f-  (5.15) 

c 


Finally,  using  fi  =  -|n|,   we  obtain 


and 


7^  - 


S  -S      CO    I  p,| 

p   ^ 

Equation  (5.I7)  Is  the  same  as  equation  (5.7).   That 

C^   ^^o 

-^  =  — 2  ^  0(5)   corresponds  to  an  equilibrium  temperature  of  the 

c    mc 

order  of  H-00   electron  volts.   At  lower  temperatures  the  effect  of 
neutral  particles  would  have  to  be  Included.   Equation  (5.I7)  would 
then  measure  the  electron  (and  hence  Ion)  density  and  not  the 
neutral  particle  density.   Equation  (5.16)  tells  us  that 
colllslonal  damping  appears  to  third  order  In  the  small  parameter 
5  . 

The  solution  corresponding  to  the  vacuum  root   s  =  -Ip   is 
obtained  by  taking  the  complex  conjugates  of  equations  (5.16)  and 
(5.17). 


C.   Plasma  Oscillations 

We  study  plasma  oscillations  in  a  region  that  has  been  covered 
extensively  in  the  literature.   In  this  region  (Ad/A)  «  1   and 

(Ad/Ac)  «  1.   The  approach  taken  here  may  be  compared  with  the 

[61 
approach  of  Lenard  and  Bernstein    who  used  a  model  most  similar 

to  ours  in  the  absence  of  a  magnetic  field  and  the  Krook  theories 

of  Lewis  and  Keller'-  •^■'  and  R.  Llboff^   .   We  show  that  the   IPP 

results  and  Libof f ' s  results  are  the  same  and  that  the  results  of 

Lenard  and  Bernstein  are  the  same  as  the  results  of  Lewis  and 
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Keller.   This  means  that  the  calculation  Is  not  sensitive  to 
collision  model  but  Is  sensitive  to  the  way  the  equatlon5are 
linearized  In  the  chosen  parameter  range.   We  shall  make  the  point 
more  explicit. 

To  make  the  computation  precise  we  choose 

r^  =  £n  «  1  (5.18) 

as   the    small   parameter.      Tlien 


^2       '       e'     '1 


=  £2  '     ^'"^  ^~^  '  ^5°^^^ 


l^ol   ^  1^0 


/\ 


I      I 


e   =£,+£_     and   e  =  e-,  +  Sp 
are  all  of  order   e-,   or  smaller  by  mass  ratios.   The  mass  ratios 
are  kept  Independent  of   £-,  .   The  remaining  parameter  q  =  Ad/A 
Is  chosen  to  satisfy 

q2  =  (Ad/A)2  =  (k2/cd2)  =  0(3')  ^  (5_20) 

The   first  hypothesis   Is   that   for     cd     >   0 

s  -   Ico   (l+h'+h'O   =   IcD   (1+a'),  (5.21) 

oxy'  o  \  -^         J 

where   h   and  h   are  supposed   0(e-,  j. 
X        y        ^-^       ^1 

The  next  hypothesis  Is  that  Hrp,  Imm:,  ^-rp     and   Irnrp^  contribute 
terms  of  order  higher  than   £-,   to  the  equation.   The  first 
hypothesis  is  verified  at  the  end  of  the  calculation,  the  second, 
which  one  verifies  by  substitution  of  (5.21)  into  the  matrix 
elements  containing  H^,  ^rnm,    H™  and   Imm^   is  not  given  here. 
With  the  two  hypothesis  the  dispersion  equation  D  =  0  may  be 
replaced  with 

-  76  - 


D 


3 


/_1  +  f^^\E^ 


\, 


i 
oy 


^1  -  ( ^  Pj 


=  0, 


(5.22) 


where 


-  ''  '       1     , 
m  ;  ^2  ^  i(l+A)  i  ' 


/  o-Ti  \ 

V 


o 


(^) 


m   2  )   m(l+A) 


^  (^2  +  iri^:AT 


g-n 


'   M   '   ,    IM 

m   ^  / 


m(l+A) 


«3=^ 


e  03 


e^(l-e-^  t)dt. 


where 


(— Jt  -  e  M.  t  +  la.  (1-e     ) 
o 


We  expand  H   In  terms  of  the  small  parameter   e   and  use 


the  Identity 


q2=^.2B'2=  (kX)  • 
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The  result  is 


H  =  1 ( 1+A )  I e 


-l(l+A)t  - 


^ot'/2  (^_^2  ^■/2)dt 


(5.23) 


With  the  substitution  y  =  q  t   and   a  =  (l+A)/q   we  get 


r 


H  =  la 


\U 


e-lay-  y  /2  ^^ 


o 


^o      2q 
^o 


(5.24) 


The  Integral  H   Is  related  to  the  Integral   I  =  [e   ^  dy 

by  use  of  the  identity 


ye' 


o  =  / -lae  o  -  


dy 


where 


y 


(5.25) 


/^  =  -lay  -  y  /2 


The  result  is 


I^  =  [ye'^Ody  =  1  -  ial^  , 


J, 


and  similarly 

00 


(5.26) 


I2  =   y^e^°dy  =  I^  -  ial^  . 
^o 

Evaluation  of   I   Is  effected  for   [crf»  1  by  means  of 

o 

the  formula  in  appendix  B  of  reference  [19].  after  a  trivial 
change  in  variables..  One  may  also  obtain  this  formula  by 
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integrating  by  parts  a  sufficient  number  of  times.   The  result 
is 


r^ 


I  = 
o 


kJ 


(5.27) 


Then 


I  o^  -  ^  (1  +  ^  +  0(\))  , 
■^     a       a  o 


ipo;  -  -^  (1  +  o(^)) 


Substitution  of  (5.28)  into  (5.24)  gives 


(5.28) 


^3  =  -  m  (l+^q?)  +  e'  - 


'o  ■ 


(5.29) 


to  lowest  order.   We  substitute  (5-29),  the  corresponding  value 
for  H^  (gotten  by  replacing  e   with  e  )   and  the  quantities 

(OzH)   (H^),  (^^),  (^^)   back  into  (5-22);  expand  to  lowest 

000        O 
order  in   ^ ' -1  5  and  equate  the  real  and  imaginary  terms  separately 

to  zero.   The  result  to  first  order  is 

2 


h  = 


3q 


X     2  ' 


y   2 


(5.30) 


or 
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s  =  la)^(  1  +  3K2/2ao^^)  -    v^/2      . 


(5.31) 


The  complex  conjugate  of  (5.31)  is  also  a  solution. 

In  the  limit  m  — ^  00  the  dispersion  equation  (5.22)  reduces 

to 

00 

2      r 

GO 


l=(v,  -^) 


o^  _s 
V-,  - 


e^dy(l-e-^iy)dy 


(5.32) 


O 


where 


^ 


s  + 


K' 


y  +  ^  (l-e-^l^) 
"1 


If  we  compute  only  to  first  order  by  the  same  techniques  leading 
to  (5.31)  we  get 


s  -  t  icD  (1  +  3K^/2co^) 


(5.33) 


with  no  damping. 

To  compare  with  Lenard  and  Bernstein    we  observe  that  if 


m 


— \     00,   if  B  =  0,   if  the  temperature  terms  are  neglected,  and 


if  one  imposes  the  contraint  that  E^    be  parallel   to  ^   then 
the  dispersion  equation  (3.23)  reduces  to  (5.32)  in  which  the 
direction  of  k  is  now  arbitrary.   The  difference  between  their 


dispersion  equation  and  (5.32)  is  that  Instead  of  the  term 

(v-,  -  0)  /s)   in  (5.32)  they  have  (-  co  /s ) .   The  missing  term  v-j^ 

stems  from  the  relaxation  effect  of  the  perturbed  velocity  on  the 
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perturbed  distribution  function  in  the  linearized  Boltzmann's 
equation  (2.13)-  The  missing  term  is  the  factor  v-,u  in  the 
quantity 

^  =  eA  ^+  v^u 


/WA 


given  in  equation  (2.13)  (in  the  limit  m  — ^  °°) '      The  term  v-,U| 
appears  when  one  linearizes  the  collision  terms  which  give  the 
macroscopic  conservation  equations.   The^"  theory  of  Lenard  and 
Bernstein  does  not  correspond  to  such  a  linearization.   The 
consequence  of  omitting  the  term  v-,u   is  that  the  root  of  the 
dispersion  formula  is  (after  one  corrects  their  equation  (38)  for 
a  typographical  error)  in  our  notation 

s  =  ±  10)^(1  +  3kV2cd^)  -  v^/2  .  (5.34) 

One  may  write  the  linearized  longitudinal  dispersion  equation 
of  Lewis  and  Keller  (equatlon/5. 7)  reference  [19])  in  our  notation 


as 


[cd^(s+v-^)+k2v^][1-(s+v^) 


e^  dy]  =  -sK^  ,  (5.35) 

^o 


where 

/  =   -(s+v^)y 1- 


Equation  (5.35)  is  the  Krook  model  version  of  the  dispersion 
equation  of  Lenard  and  Bernstein.   The  result  of  computing  the 
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root  corresponding  to  the  plasma  oscillation  is  precisely  equation 

(5.3^)-   This  reflects  the  fact  that  the  linearized  Boltzmann's 

equation  of  Lewis  and  Keller  lacks  the  same  term  as  does  the 

equation  of  Lenard  and  Bernstein.   If   one  computes  from  Libof f ' s 

equations    ,  on  the  other  hand,  one  should  obtain  our  equation 

(503)'   This  reflects  the  fact  that  Libof f ' s  equations  mimic  the 

Boltzmann's  equations  and  when  linearized  contain  a  term  like   v-,u  . 

J.-— I 

In  the  chosen  parameter  region  the  mechanism  of  damping  to 
first  order  is  given  by  the  ion-electron  collision  frequency   Vp   as 
seen  in  equation  (5.3l)-   A  properly  linearized  Krook  theory  should 
give  the  same  answer  in  the  same  parameter  region  when  the  collision 
frequencies  are  suitably  identified. 

Equation  (5.31)  agrees  with  the  results  of  Johnston  and 

[211 
Shkarofsky     to  the  order  given.   Their  results  are  computed  from 

a  set  of  linearized  macroscopic  equations  obtained  by  estimating 

collisional  effects  from  the  original  Fokker-Planck  coefficients  to 

lowest  order. 

Their  equations  also  satisfy  conservation  of  mass,  momentum 

and  energy  for  finite  ion  mass.   It  is  therefore  not  surprising 

that  to  the  order  of  equation  (5-3l)  the  results  agree. 
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6.   Summary  and  Conclusions 

The  actual  Boltzmann's  and  Fokker-Planck  equations  are  too 
complicated  at  present  to  be  of  use  In  constructing  a  dispersion 
equation.   For  this  reason  one  seeks  collision  models  which  retain 
as  many  of  the  properties  of  the  original  equations  as  possible. 
There  are  tvjo  models  which  retain  the  property  of  giving  the 

macroscopic  conservation  equations:  Llbof f ' s  generalization  of  the 

[pl 
Krook  model    and  the  isotropic  Pokker-Planck  model  Introduced 

here.   Of  the  two  we  expect  that  the  IFF  model  is  superior  since  it 

retains  the  physical  property  of  describing  a  diffusion  process  in 

velocity  space.   ¥e  have,  therefore.  Increased  confidence  in  those 

results  which  are  the  same  in  both  theories  and  where  there  Is 

disagreement  we  expect  the  IFF  results  to  be  better. 

In  general  one  expects  that  long  wavelength  disturbances 
behave  macroscopically.   For  these  disturbances,  the  magnetohydro- 
dynamlc  and  plasma  electromagnetic  modes,  there  is  little  difference 
in  the  models.   Differences  appear  as  small  higher  order  .       .    ■ 
corrections.   In  section  5  where  the  Faraday  rotation  and  plasma 
oscillations  are  studied  the  two  theories  give  the  same  corrections 
since  the  collision  frequencies  are  chosen  to  be  small  with 
respect  to  the  representative  frequencies. 

The  greatest  disparity  between  the  Krook  and  IFF  models 
appears  in  the  results  for  the  microscopic  Larmor  resonances  in 
the  long  wavelength  limit.   We  see  in  sections  4B  and  4D  that  the 
Krook  model  does  not  describe  a  diffusion  process  in  velocity 
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space  since  It  gives  essentially  the  same  results  as  the 
collisionless  model.   Where  the  Krook  model  predicts  an   infinite 
number  of  Larmor  resonance  modes  only  at  propagation  perpendicular 
to  the  applied  magnetic  field  the  IPP  model  predicts  an  infinite 
number  of  Larmor  resonances  at  arbitrary  directions  of  propagation 
(except  parallel  to  the  applied  field).   At  perpendicular 
propagation,  for  each  value  of  the  integer  N  multiplying  the 
Larmor  frequency  there  is  but  one  Krook  resonance  and  an  infinite 
sequence  of  IPP  resonances  differing  by  integral  multiples  of  the 
collision  frequency.   This  last  difference  reflects  the  mathematical 
property  that  the  IPP  differential  operator  has  an  unbounded 
discrete  spectrum.   The  unboundedness  of  the  spectrum  is  reflected 
again  at  parallel  propagation  where  there  are  at  the  Larmor 
frequency   (|n|  =  l)   an  infinite  number  of  resonances  which  differ 
by  integral  multiples  of  the  collision  frequency.   Such  modes  are 
not  found  from  the  Krook  theory.   Here  the  mean  free  path  was  small 
with  respect  to  the  wavelength.   Por  mean  free  path  large  with 
respect  to  wavelength  we  find  in  section  4F  strongly  damped 
resonances  at  the  Larmor  frequency  again  at  parallel  propagation. 

In  section  hC   we  study  the  new  Larmor  resonance  modes  at 
arbitrary  angles  of  propagation  in  the  non-uniform  collisionless 
limit  obtained  by  keeping  the  ratio  of  mean  free  path  to  wavelength 
fixed  and  small  as  the  collision  frequencies  become  small.   We 
find  that  the  new  Larmor  resonance  modes  remain  in  this  limit. 
These  results  indicate  the  possibility  of  experimental  verification 
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if  one  viewed  radiation  from  a  plasma  at  angles  other  than 

perpendicular  to  the  applied  magnetic  field.   They  are,  however, 

at  present  only  suggestive  since  we  have  not  examined  the  mechanism 

for  exciting  such  modes. 

[221 
S.  Tanaka  et  al.     have  studied  microwave  radiations  from 

a  dc   discharge  plasma  In  a  magnetic  field.   Although  their 

measurements  appear  to  be  of  Larmor  resonances  their  results  are 

not  sensitive  enough  to  distinguish  between  models, especially  since 

they  view  the  radiation  at  a  direction  perpendicular  to  the  applied 

magnetic  field.   Their  measurements  do  agree  with  the  rough 

observation  that  for  a  denser  plasma  the  Larmor  resonances  are 

given  by   s  =  tN|n[  +  •••   where   |n[  >  1   and  N  =  t  1   Is  not 

a  root  of  the  dispersion  equation.   As  the  density  becomes  small 

the  roots   s  =  t   i[^I   are  approached.   This  observation  Is 

fPl 

contained  In  Llboff 's  study   ■"  since  those  roots  which  approach 

s  =  _  llfil   as  CO   and   v  both  become  small  linearly  with  the 
density  are  among  his  lowest  order  plasma  electromagnetic  modes. 
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Appendix  A 

List  of  Integrals: 
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Q-^   -    e 
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Q    =  g-vy  Li(^+%)(i_e-(v+ifi))2  ^  g-i(f+fiy)^_ 
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